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Introduction 


Introduction 


This handbook is a reference that you can take into the MST224 exam, 
and it may be valuable when you start to apply your knowledge in other 
modules. It will be more effective if you are already familiar with it before 
you sit the exam, and we suggest that you consult it when you attempt the 
assignments in the module. This handbook is not designed as a teaching 
document, and reading it is not a substitute for studying the module units. 


The first few sections consist of general mathematical reference material, 
largely based on the topics in Unit 1. It is not intended to be a 
self-contained or logically complete account of basic mathematics; it is just 
a set of definitions and results. Some of these results are used repeatedly in 
the module. You will also find some additional results or definitions that 
are not covered in Unit 1, but are useful to have in a mathematical 
reference booklet. 


The later sections of the handbook are brief summaries of the units, 
emphasising the most important results. 


1 Notation 


1.1 Greek alphabet 

a A alpha e I iota p P rho 

B B beta kK K_ kappa o D sigma 
y TI gamma A A lambda 7 T tau 

6 A delta uw M mu v YT upsilon 
e E epsilon vy N ou ® phi 

¢ Z zeta €é 2B x x X chi 

n H eta o QO. omicron wv W psi 

6 O theta a Il pi w Q omega 
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1.2 Symbols 


is equal to 

is not equal to 

is approximately equal to 
plus or minus 

minus or plus 

less than 

less than or equal to 
greater than 

greater than or equal to 
positive square root 

the number 2.718 28... 
the number 3.14159... 
infinity 

the integers 

the real numbers 

the complex numbers 


QAN BA SIV VIAAH HR I 


1.3 Limits and sums 


An ordered list of numbers 29,71, #2, .--, Lp,... is said to converge to the 
limit x if successive terms in the ordered list are better and better 
approximations to x. We write ‘7, > x as n > ov’ or ‘Tim Ba =o’. 
Given numbers aj,a2,..., Qn, we define 
n 
Soa to mean a,+a2+---+@n-1+4n. 


i=1 


2 Numbers 


2.1 Real numbers 


The integers are the positive and negative whole numbers, together with 
zero. Non-integer numbers that can be expressed exactly as fractions are 
called rational numbers; those that cannot be so expressed, such as V2, 
e and 7, are called irrational. The collection of the rational numbers 
(including the integers) and the irrational numbers is called the set of real 
numbers. 


When a real number is expressed in decimal notation, if it is approximated 
then the approximation can be given to so many decimal places, or to so 
many significant figures. For example, 1.4142 is the approximation to 
V2 to four decimal places and five significant figures, while 0.000000 342 is 
given to nine decimal places but three significant figures, and 342000 has 
no decimal places but at least three significant figures. The process of 


reducing the number of decimal places or significant figures to which a 
number is expressed is referred to as rounding. To round a given number 
to n decimal places or n significant figures, take the number expressed to 
n decimal places or n significant figures that is closest to the given 
number, where it is conventional to round the digit 5 up. (Be aware that 
other conventions for rounding exist, and that computer programs and 
calculators do not always use this convention.) 


A number given in the form +b x 10°, where 1 < b < 10 and c is an 
integer, is said to be in scientific notation. For example, the number 
342000 can be expressed as 3.42 x 10° in scientific notation. 


For any two real numbers a and b with a < b, we write: 
{a, | for the set of all real numbers x such that a < x < b: 
{a, b) for the set of all real numbers x such that a < x < b; 
(a, | for the set of all real numbers x such that a < x < b; 
(a,b) for the set of all real numbers x such that a < x < b. 


These sets of numbers are called intervals. The interval [a, }] is a closed 
interval, (a,b) is an open interval, and the two others are half-open 
intervals. For practical problems the distinction between open and closed 
intervals is rarely significant. If we know that a real number z is 1.274 to 
three decimal places, then x lies in the interval [1.2735, 1.2745). 


2.2 Complex numbers 


A complex number 2 is written in Cartesian form as z = a + bi, where 
a and bare real numbers, and i? = —1. We refer to a as the real part 
of z, written Re(z), and to b as the imaginary part of z, written Im(z). 


Complex numbers can be added, e.g. 
(a+ bi) + (c+ di) = (a+c) +(b+d)i, 
or multiplied, e.g. 
(a + bi)(c + di) = (ac — bd) + (ad + be)i. 
These formulas make use of the ordinary rules of algebra, together with 
the relation i? = —1. 
The complex conjugate of z = a + bi is 7 =a — bi. Note that 
22 = (a+ bi)(a— bi) =a? + 0? 
is a positive real number (unless a = b = 0). The modulus of z is the 


number |z| = /zz = Va? +B. 


To calculate a quotient of complex numbers, multiply top and bottom by 
the complex conjugate of the bottom, e.g. 


at+bi _(at+bi)(c—di) _ (act+bd be — ad i 
c+di  (c+di)(c—di) C+é C+) 
which expresses the quotient in the form p+ gi with p and q real. 


2 Numbers 


The notation 2* is also used in 
some texts for the complex 
conjugate of z. 


= and 2° are read as ‘z bar’ and 
“z star’. 
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Polar coordinates are discussed 
in Subsection 5.2. 


The Argand diagram is a representation of complex numbers as points 
in a plane, where the complex number a + bi is represented by the point 
with Cartesian coordinates (a,b). A point can also be represented in polar 
coordinates as (r,@) and related to its Cartesian coordinates by 


a=rcos#, b=rsin@. 


The angle coordinate @ is referred to as an argument of z, written arg(z), 
and the unique value of @ in the range —7z < 6 < 77 is referred to as the 
principal value of the argument, written Arg(z). 
The multiplication of complex numbers in polar form is given by the rule 
(r,0) x (s,¢) = (rs,@+ 4). So powers of a complex number can be 
expressed as (r,@)” = (r",n8), for n a positive integer. The special case 
when r = 1 is known as de Moivre’s theorem: 
(cos @ + isin 0)” = cos(n@) + isin(n8). 
Euler’s formula is 
e® = cos6 + isin#. 
This extends to any complex number a + bi as 
ettbi — eth — 6*(cosh + isind). 
The exponential form of a complex number z = r(cos6 + isin @) is 


z=re™. 


This form is useful for multiplying, dividing and taking powers of complex 
numbers. 


3 Functions and graphs 


3.1 Functions 


A variable is a quantity, represented by a symbol, that can vary over a set 
of values. If its value does not vary, then it is a constant. 


Any expression or formula that involves a variable x, and whose value is 
uniquely determined by the value of z, is called a function of wx. 


If a variable y is a function of z (i.e. if y is equal to a function of x), then 
we call x the independent variable and y the dependent variable, and 
we may write y = y(x). Here y(a) stands for the function of x (i.e. for the 
formula involving x). 


If f and g are two functions, then their sum is a function f + g defined by 
(f + 9)() = f(x) +9(a) (for all x). 


Moreover, if A and B are any two numbers, then the function Af + Bg is 
defined by 


(Af + Bg)(x) = A f(x) + Bg(x) (for all x). 


3 Functions and graphs 


The function f(g(x)) is called the composite function or composition 
of the functions f and g. 


The graph of a function f(x) is the curve in the ry-plane whose equation 
is y = f(z). 

A continuous function is one whose graph has no breaks or jumps in it, 
i.e. it can be drawn without lifting your pen from the paper. 


A constant function f(x) is one that assigns the same value to any 
input. Its graph is a straight line parallel to the a-axis. A special case is 
the zero function, which assigns the value 0 to any input. 


A linear function is one having the form a,x + ag (with a; # 0), where 
a, and ao are constants. Its graph is a straight line with slope a, and 
y-intercept ap. 


A quadratic function is one having the form a2x? + ax + a9 (with 

az # 0), where a2, a), a9 are constants. Its graph is a parabola, similar in 
shape to the one shown in Subsection 5.3 if ag > 0, but the other way up if 
ag <0. 


A cubic function is one having the form a3z° + agx? + a,x + ap (with 
a3 # 0), where a3, 2,4 ,a9 are constants. 


3.2 Polynomials 


Linear, quadratic and cubic functions are all particular examples of 
polynomial functions, or simply polynomials. An nth-order polynomial, 
or polynomial of degree n, is a function of the form 


nx" + an12") +--+ +012 +49, 


where n is a positive integer, x is a variable and ao, ai,..-, Qn, are 
constants with a, #0. A linear polynomial has n = 1 (orn =0),a 
quadratic polynomial has n = 2, and a cubic polynomial has n = 3. 


The roots of a polynomial p(x) are the solutions of the equation p(x) = 0. 
Every polynomial of degree n can be written as a product of a, and 

n factors of the form a — cy (k = 1,2,...,n), with each c;, a complex 
number (which may be real). Each of these factors corresponds to a root 

a = cy of the polynomial. If a factor « — ¢ occurs more than once, then the 
root x = c is a repeated root; repeated roots are also sometimes referred 
to as equal roots or coincident roots. 


The roots of a quadratic equation ax? + bx + ¢ = 0, a 4 0, are given by the 
formula method as 


—b + Vb? —4ac 
2a . 


The quantity 6? — dac is referred to as the discriminant of the quadratic 
equation. 
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Some texts refer to y = ba” as 
an exponential function; this is 
not to be confused with the 
exponential function, exp x. 
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To factorise a polynomial is to express it as a product of two or more 
polynomials of lower degree. For example, the difference of two squares 
x? — a? factorises as 


x” —a? =(x—a)(x+a), 


and the perfect square x? + 2axr + a? factorises as 


x’ + 2ax +a? = (2+a)(x+a) = (4 +a). 


3.3 Exponentials and logarithms 


A function of the form y = ba*, where a and 6 are non-zero constants (with 
a> 0 and a #1), and z is real, is said to exhibit exponential 
behaviour. In a’, a is referred to as the base and x as the exponent (or 
index or power). Properties of such a function include 


a*=1/a", 

a” x a =a", 

a®/a¥ =a", 

(a*)" = a" = (a")*. 
The function e*, where e = 2.718 28..., is referred to as the exponential 
function. It is also written as expz. 


The natural logarithm function In x is defined to be the inverse 
function of the exponential function exp z, i.e. each reverses the effect of 
the other, so that 
In(expx) =a for all real x, 
exp(Inz) = for all real x > 0. 
In other words, if e” = x, then y = Ina, and vice versa. 
The natural logarithm function In x, for x > 0, has the properties 
Inl=0, 
In(1/x) = —Inz, 
In(ay) =Inz+Iny, y>0, 
In(a/y)=Inx—Iny, y>0, 
Ina’ =ylng. 
Any function y = ba” can be written in the form y = be**, where k = Ina. 


Another logarithm function is logy x, for > 0, where y = logy x if 
10” = x (and vice versa). The properties given above for In also hold for 
logio- 


Although they are not used in this module, we note for general reference 
purposes that the hyperbolic functions sinh x, cosh and tanh x are 
defined as combinations of exponential functions: 


sinha = 3(e* —e”*), 


cosh = $(e7 +e), 


tanh zc = —— = ——_.. 
coshaz e™ +e-* 


The inverses of these functions are 


arcsinh x = In(x + V2? +1), 
arccoshx = In(r+ V2?-1), x>1, 


1 
arctanh x = +in( *=) ¢ lel<k 
a 1-z 


3.4 Graphs of some common functions 


3 Functions and graphs 


yy yy ys 
Ea a3 Vr 
> > > 
zr a x 
yt yy 
1 Ina, 
x 
> > > 
x x x 
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4 Trigonometry 


4.1 Radians and degrees 


In this module we usually measure angles in radians rather than degrees. 
There are 27 radians in a full circle, corresponding to 360°, so 1 radian is 
(180/7)° ~ 57°. An advantage of working in radians is the simplicity of the 
formula for the are length subtended by an angle in a circle of radius r: 
the length of the arc subtended by an angle of @ radians is simply r@. The 
following radian measures of standard angles are worth knowing: 


e aright angle is } radians 
e the angles of an equilateral triangle are § radians. 


An angle is acute if its radian measure lies between 0 and oa and obtuse 
if its radian measure lies between 5 and 7. 


4.2 Trigonometric functions and their inverses 


For an acute angle 0, the values of trigonometric functions are related to 
the ratios of lengths of the sides of a triangle as follows. 


Function Definition for acute angles Definition in terms 


in terms of triangle shown of sin and cos 
sin? £ 
G 
b 
cos 6 = 
c 
a sin@ 
tand = 
st b cos 
cot @ 8 a Y 
a sin@ 
c 1 
tf) = 
= b cos @ 
1 
cosec # g a 
a sin# 


4 Trigonometry 
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Inverse function Definition 


arcsin © = 6 where sin@ = x 
arccos © = 6 where cos# = x 
arctan © = 6 where tan@ = x 
arccot © = 6 where cot? = x 
arcsec x = 0 where sec@ = x 
arccosec © = 6 where cosec @ = x 


4.3 Two useful triangles 


From the two triangles in the margin it can be seen that: 


ine = 2 zav3 pe 
sing = 3: cosg =p: tan§ = 
eas x 2 Hi fh: 
cosec § = 2, sec § = Ys: cot = V3; 
i 1 J a 
sint =z, cost = 7B tanZ=1, 
cosec | = V2, sect = V2, cot $=1; 
sing = 8, cos 3 = $, tan ¥ = v3, 
rai sec Z = zt 
cosec 3 = 7, sec § = 2. cot z= 7 
Other values of the trigonometric functions worth remembering are: 
ig! ng 
sin0 = 0, cos0 = 1, tan0 = 0; 
x x0. 
sn 5 =1 cos 5 = 0; 
sinz = 0, cos7 = —1 tan z= 0; 
sin se =-l, cos 32 = 


4.4 Trigonometric identities 


Pythagoras’s theorem states that for any right-angled triangle, if c is 
the length of the hypotenuse (the side opposite the right angle) and a and 
b are the lengths of the other two sides, then 


e=84B. 
This leads to the following trigonometric identities: 
sin? 6 + cos” @ = 1, 
tan” @ + 1=sec? 0, 
1+ cot? 0 = cosec? 0. 


Addition formulas 
sin(a + 8) = sinacos § + cos asin f, 
sin(a — 8) = sinacos 8 — cosasin B, 
cos(a + 8) = cos acos § — sina sin 3, 
cos(a@ — 8) = cos acos 8 + sina sin 8, 

tana +tanf 

1—tanatan §’ 

tana — tan 

T+ tanatan 


tan(a + 8) = 
tan(a ~ 8) = 


sinacos § = 5sin(a+ 3) + $sin(a — 8), 


3 
cos asin 8 = $sin(a + 8) — $sin(a— 8), 
cos acos B = $ cos(a +£B)+ $ cos(a — 6), 
sinasin 8 = $.cos(a — 8) — $ cos(a + 8). 
In particular, these formulas give 
sin(a +27) =sina, cos(a+2m)=cosa, tan(a+7) = tana; 
sin(—a) =—sina, cos(—a) = cosa, tan(—a) = — tana. 


Double-angle formulas 


sin 2a = 2sina cosa, 


ae 2 
cos 2a = cos? a — sin? @ = 1 — 2sin? a = 2cos* a — 1, 


ace == 2tana 
~ 1=tan? a’ 
sin? a = $(1 —cos 2a), 


cos? a = 31 + cos 2a). 


Cosines of related angles 
cos(§ —a)=sina, cos($ +a) =—sina, 
cos(7 — a) =—cosa, cos(7+a) = —cosa. 


4.5 General sinusoidal functions 
A sinusoidal function or sinusoid is a function x(t) of time ¢ of the form 
x = 29 + Acos(wt + ¢) = x9 + Asin(wt + ¥), 


where 2p is a constant, A is a positive constant called the amplitude, w is 
a positive constant called the angular frequency, and ¢ and v are 
constants called phase constants. 


4 Trigonometry 
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A sinusoidal function oscillates between xp — A and xp + A, repeating the 
same pattern of oscillations through each time interval of length 27/w, 
known as the period of the function. For these reasons, sinusoidal 
functions are examples of oscillatory functions and of periodic functions. 


The phase constants in the two forms of the sinusoidal function are related 
according to ¢ = = + (or equivalently ¢ = y — 3). 


Alternative forms of the sinusoidal functions are given by 


x = 29 + Acos(wt + ¢) = x9 + Bcos(wt) + Csin(wt), (1) 
& = 29 + Asin(wt +) = x9 + Dsin(wt) + Ecos(wt). (2) 


In equation (1) we have B = Acos@ and C = —Asin@. Solving these 
equations gives A = VB? + C? and ¢ = arccos(B/A). 


In equation (2) we have D = Acos~ and E = Asin. Solving these 
equations gives A = VD? + E? and w = arccos(D/A). 


5 Geometry 


5.1 Cartesian coordinates 


The Cartesian coordinates (, y) of a point P in a plane specify the 
position of that point relative to two perpendicular axes, the z-axis (or 
horizontal axis) and y-axis (or vertical axis), which meet at a point O 
called the origin, with Cartesian coordinates (0,0). The directions of the 
axes indicate increasing numerical values for the x- and y-coordinates. 
Values of x to the right of the y-axis are positive, and those to the left are 
negative; similarly, values of y above the z-axis are positive, and those 
below are negative. The four parts into which a plane is divided by 
Cartesian coordinate axes are known as quadrants of the plane. A plane 
on which Cartesian coordinate axes have been specified is often referred to 
as the zy-plane. 


y 
y 
P Second quadrant | First quadrant 
ae o(x,y) 2<0 x>0 
y>0 y>o0 
> 
ee Oo 2 
origin r<0 «>0 
O = y<0 y<0 
x-axis Third quadrant | Fourth quadrant 
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5.2 Polar coordinates 


The point P whose polar coordinates are (r,@) has Cartesian coordinates 
(x,y) where 


x=rcos#, y=rsiné. 


The value of r is always positive (except at the origin, where it is zero). 
For a given point P, the value of @ is not unique: we can add or subtract 
any integer multiple of 27 and obtain another value for @ that describes the 
same point. The value of @ satisfying —7 < 0 < 7 is called the principal 
value of 6. 


The relationship between 
polar and Cartesian 
coordinates 


5.3 Plane figures and curves 


A closed plane figure with straight sides is called a polygon. A polygon 
with 3 sides is a triangle, one with 4 sides is a quadrilateral, one with 
5 sides is a pentagon, one with 6 sides is a hexagon, and in general one 
with n sides is called an n-gon. 


A polygon is said to be regular if all its sides have equal length and all its 
angles are equal. A regular triangle is referred to as an equilateral 
triangle, and a regular quadrilateral is a square. 


An isosceles triangle is one with two sides of equal length (or 
equivalently with two equal angles). A right-angled triangle is one in 
which one angle is a right angle. The angle sum of a triangle is 

m radians (180°). 


A parallelogram is a quadrilateral with opposite sides parallel. 

A rectangle is a parallelogram all of whose angles are right angles. 

A square is a rectangle all of whose sides have equal length. The angle 
sum of a quadrilateral is 27 radians (360°). 


The area of a triangle = $ x base x height. 


The area of a parallelogram = base x height. 


height height 


——— 
base base 


The areas of a triangle and parallelogram 


iz 
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A circle is a set of points in a plane that are a constant distance from a 
fixed point in the plane. The fixed point is the centre of the circle, and 
the constant distance is its radius. If a straight line cuts a circle at two 
points, then the segment of that straight line within the circle is known as 
a chord of the circle. The length of a chord that passes through the centre 
of a circle is the diameter of the circle. The terms diameter and radius 
are also used to refer to a chord through the centre of a circle and to a 
straight line from a point on the circle to its centre, respectively. 


angle subtended 
by arc 


Each continuous segment of a circle is known as an arc of the circle; the 
angle made at the centre of a circle by two radii drawn from the ends of the 
arc is known as the angle subtended by the arc. (In such circumstances 
the arc itself is sometimes called the arc subtended by the angle.) 


The distance around a circle is known as its circumference, and for a 
circle of radius r this is given by 27r. The area of a circle of radius r 
is xr?. The arc length of an arc of a circle of radius r subtended by an 
angle @ radians is r@. The area of a sector of a circle of radius r 
subtended by an angle @ radians is $170. 


The equation of a circle in the ry-plane with centre (a,b) and radius r 
is (x — a)? + (y—6)? =r?. Other curves in the ry-plane that can be 
represented by quadratic formulas are the ellipse, the parabola and the 
hyperbola, examples of which are shown below. 


nat Nh Pt 


Ellipse Parabola Hyperbola 


A straight line between two distinct points on a curve is known as a chord 
of the curve. A straight line that just touches a curve is known as a 
tangent to that curve at the point where it touches. 


6 Differentiation 
6 Differentiation 


6.1 Notation and terminology 


If f(z) is a function, then its derived function or derivative f'(x) is 
defined by 


f(e+h)— f(e), 


a eee: 


The process of calculating f’(x) from f(z) is called differentiation of 
f(x) with respect to x. Differentiation with respect to x can also be 


d 
denoted by the symbol E written to the left of the expression or variable 
le 


being differentiated, so that, for example, 


d (cx) . 
af) and or both mean f’(z). 
If y = f(x), then f’(x) can also be written as y/ or oe where to save space 
d. 
we often print dy/dx in place of = 


When the independent variable is t (time), we often use a dot to indicate a 
derivative, so that % means the same thing as u’(t) or du/dt. If x is the 
position of a particle along an axis, then dx/dt or « represents the 
velocity along the axis, and d*x/dt? or # represents the acceleration 
along the axis. 


The notation f’(x) is referred to as function notation, dy/dx as Leibniz 
notation, and % as Newton’s or Newtonian notation. 


The derivative of a derivative is called a second derivative. For example, 
the second derivative of the function f(x), denoted by f”(x), is the 
derivative of f’(x). If y= f(x), then the second derivative is also written 
as y” or d?y/dx?. If u is a function of t, then its second derivative can be 
written as ti. 

Third and higher derivatives are defined analogously. The nth 
derivative of f is denoted by f or, if y = f(x), by d"y/dx"; n is referred 
to as the order of the derivative. The prime and dot notations are not 
ed fe higher derivatives, except that f’” is sometimes used in place 

of fF, 


A complex-valued function f(x) = g(a) +ih(«), where g and h are 
real-valued functions, can be differentiated in a natural way as 


f'(x) =g'(z) +ih'(z). 
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6.2 Rules of differentiation 


Constant multiple rule: If k is a constant and u is a function of x, then 


d du 
ee = ipod i —(ku) = k— 
(ku)’ = ku’, or equivalently an (ke) k in 
Sum rule: If u and v are functions of x, then 
d du. dv 
sgl Patt ivalently = +o. 
(u+v)!=u'+v', or equivalently (u+v) ae + ae 


Product and quotient rules: If u and v are functions of x, then 


du _ dv 


(uv)! =u'v+uv', or equivalently  (w) = 


and 


v 


or equivalently <(*) = (3 a +) Je. 


This is sometimes called the Composite rule or chain rule: If g and wu are two functions, and 
‘function of a function’ rule. h(x) = g(u(x)), then 


h(a) = g/(u(x)) w'(z). 
Another way to write this is as 

dh _ dh du 

dx du dx" 
Implicit differentiation: Given an equation connecting two variables x 
and y, we can use implicit differentiation to calculate dy/dx by 
differentiating both sides with respect to x and then solving the resulting 


equation algebraically for dy/dx (instead of solving for y before 
differentiating). For example, 


(1312) — 9724? + 9x3, Y 
gt y) ss +2n"y 7: 
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6.3 Standard derivatives 


In each case, a is a constant. 


Function Derivative 
a 0 
a ar’! 
oe ae“ 
In(azx) x 
x 
sin(ax) acos(ax) 
cos(ax) —asin(ax) 
tan(az) asec?(ax) 
cot (ax) —acosec*(ax) 
sec(aa) asec(ax) tan(ax) 
cosec(ax) —acosec(ax) cot(ax) 
s a 
arcsin(az) Vitae 
a 
arccos(ax) Fras 
a 
arctan(ax) Dee 
arccot(az) oe 
1+ a?x? 
P a 
arcsec(aa’) jae 
a 
arccosec(ax) Teale Veni 


6 Differentiation 
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6.4 Stationary points 


The gradient of a function f at a point zo is the slope of the tangent to 

the graph of f at that point, and is given by the derivative of f at that 

point, i.e. f’(zo). A function is increasing on an interval if its gradient is 

positive throughout that interval; it is decreasing if its gradient is 

negative throughout that interval. A stationary point of f is a point x 

where the gradient is zero, i.e. f’(x9) = 0. 

A function is smooth if it is continuous and has a continuous derivative. 

Consider a stationary point x9 of a smooth function f. 

e zo is a local maximum if, for all x in the immediate vicinity of xo, 
f'(x) > O if  < x and f(x) < 0 if x > zo. An alternative condition is 
f" (x0) < 0. 

e = xo isa local minimum if, for all x in the immediate vicinity of xo, 
f'(x) < Oif x < x and f'(x) > 0 if z > zp. An alternative condition is 
f" (ao) > 0. 

e zo is a point of inflection if f” (xo) = 0 and f”(x) changes sign as x 
increases through zo. 

A global maximum of a function f is a point xo such that f(x) > f(x) 

for all 2 where f is defined. A global minimum of a function f is a 

point x9 such that f(29) < f(x) for all x where f is defined. A function f 

is bounded above by an upper bound A if f(x) < A for all x where f 

is defined. A function f is bounded below by a lower bound B if 

f(x) > B for all x where f is defined. 


6.5 Curve sketching 
The following is a possible procedure for sketching the graph of y = f(z), 
where f(x) is some given function. 


1. Check whether f(z) is a standard function whose graph you already 
know, or is a simple modification of such a function. If not, proceed to 
Step 2. 


2. Determine how y behaves when = is very large and positive, and when 
zx is very large and negative. 


Look for any obvious symmetries or repetitions in the behaviour of f. 
Find where the curve crosses the x- and y-axes, if at all. 


Look for any values of 2 at which f(a) is undefined, and examine the 
behaviour of f(x) near these values of x. 

6. Find the locations of any local maxima, local minima or points of 
inflection. 

7. Try to determine whether there are any intervals over which the 
function is increasing or decreasing. 


6 Differentiation 


8. Transfer the information found in Steps 4 and 6 to a sketch graph, 
then use this information together with any information found in 
Steps 2, 3, 5 and 7 to try to sketch a smooth curve. If you are still 
unsure about any parts of the curve, choose suitable values of « and 
plot the corresponding points (x, f(x)) before completing the curve. 


Lines y = c where f(x) — c as x — +00, and lines x = c where 
f(x) + +00 as x +, where c is a constant, are known as asymptotes of 
the graph of the function f(x). 


6.6 Taylor polynomials and series 


Factorials 
For any positive integer n, we define n factorial, written n!, by 
nl! =1x2x3x---x(n—-1)xn. 


The first few factorials are 1! = 1, 2! = 2, 3! = 6, 4! = 24. We also define 
ol=1. 


Taylor polynomials 


For a function f(x) with n continuous derivatives near « = a, the Taylor 
polynomial of degree n about x = a or the nth-order Taylor 
polynomial about x = a is 


pal) = f(a) + («~a)f"(a) + (0 — a? f"(a) +--+ Glo a)" F(a). 


When used to approximate f(x) near x = a, we refer to this polynomial as 
the nth-order Taylor approximation to f(x) near x =a, and write 


fle) = f(a) +(e a) f(a) + le ~ a)? f"(a) +--+ Sea)" Fa), 
In particular, n = 1 gives the tangent approximation 
J (x) = f(a) +(e —a)f"(a), 
and n = 2 gives the quadratic approximation 
a) & f(a) + (a —a)f'(a) + 5(« — a)? f"(a). 
These approximations are good when z is close to a. 
Taylor series 


The Taylor series about x =a for a function f(x) with infinitely many 
continuous derivatives near x = a is 


fle) = f(a) +(e a) f(a) + (2 —a)?f"(a) + 


A Alaa)" f(a) shay 
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Some standard Taylor series about « = 0 


Sees Us ee na1__1 2n-1 
sing =a2— 2% + ae —---+(-1) @n—1)" feet, 
gt 12 lu n_1_oon 
costs = 1— Fr +z —-++++(-1) Bi feisty 
x 1» ln 
P=H1l+e+sat+---+—a tes, 
2! n! 
121s n-il in 

In(l+2)=2—-5¢ +3 = ++e(-1) aor (-l<a<l). 


Small-angle approximations 


If the angle @ is small (sometimes written << 1) and is measured in 
radians, then we can obtain useful approximations by truncating the above 
Taylor series: 


sin ~ 0, 
cosé ~ 1— 467. 


7 Integration 


7.1 Notation and terminology 


The indefinite integral of a continuous function f(x) is 


[tea = F(a) +C, 


where F is a function such that F’(x) = f(x), known as an integral or 
antiderivative of f, and C is a constant, often referred to as an 
arbitrary constant or constant of integration. 


The definite integral of a continuous function f(x) from a to b is 


b 
| J (x) dx = [F(x)]’ = F(b) — F(a), 


where F is any integral of f. The numbers a and b are called the lower 
limit of integration and upper limit of integration, respectively. If 
the areas bounded by the graph of f(x) above and below the z-axis 
between a and b are A; and Ag, respectively, then 


[souna — Ao. 


7 Integration 


The process of finding an indefinite or definite integral is known as 
integration, and the function f being integrated is known as the 
integrand. If F; and F» are two integrals of f, then they differ by a 
constant, i.e. Fy (a) = F(x) + C, where C is a constant. 


7.2 Rules of integration 


Constant multiple rule: 

[rso) dx = k [ #2) dx (where k is a constant). 
Sum rule: 

[6 +06) ar = f sear + f o(a)ae. 
Integration by substitution: 

/ (ola)) (a) dz = { flu)du (where w= 42), 
or in Leibniz notation 


[eo Gae= f su au. 


The following formula, which can be derived by integration by 
substitution, is also useful: 


[2 dx =I\n|g(x)|+C (where g(x) 4 0). 


Integration by parts: 
[ e)a@ae= se) a2) f 1(@) a(a)ae. 
For definite integrals, 


b b b 
[ sos @ae= (ge) ae))t = [ f@)afe) a. 


A function f is odd if f(—r) = —f(x). If f is odd and a is a positive 
constant, then 


‘ f(x) dx =0. 


7.3 Standard integrals 


In each case, a is a non-zero constant, b is any constant, and n is any 
integer. When using the table below to obtain indefinite integrals, add an 
arbitrary constant. 
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Function Integral 
a ar 
got 
‘i -1 
o* (a#—1) a+l1 
1 pl 
-hh ti) 
ax+b a lent) 
- So 
a 
In(az) x(In(ax) — 1) 
1 
sin(az) > cos(ax) 
de 
cos(ax) —sin(axr) 
a 
1 
tan(ax) —7 In|eos(az)| 
Wut 
cot(ax) a In |sin(az)| 
1 
sec(ax) - In |sec(ax) + tan(az)| 
cosec(axr) = In |cosec(ax) — cot(azx)| 
a 
2 1 
sec? (ax) i tan(axr) 
2 1 
cosec*(ax) =a cot(ax) 
; Mas 
xsin(ax) 3 (sin(ax) — ax cos(ax)) 
1 
xcos(axr) z (cos(ax) + ax sin(ax)) 
1 1 © 
a? +a? a sass @) 
1 1 a= % 
at 
(a —a)(x —b) a—b |x-—b 
1 
[>= In(x + Vx? + a? 
V2? +a? ( ) 
1 
= In|x + Vx? — a? | 
‘x? — a? 
1 = (2 
ao arcsin (=) 


Unit 2 summary 


Unit summaries 


Unit 1 Getting started 


The material in Unit 1 is summarised in Sections 1-7 of this Handbook. 


Unit 2 First-order differential equations 


ds 


A differential equation is an equation that relates an independent 
variable x, a dependent variable y, and one or more derivatives of y. 
Its order is the order of the highest derivative that appears. So a 
first-order differential equation for y = y(x) contains its first 
derivative, dy/dx, but no higher derivative of y. This unit considers 
only first-order differential equations that can be expressed in the form 


dy 
dp FY) 


where f(2,y) stands for some expression that may contain either or 
both of the variables x and y. 

A solution of a differential equation is a function y = y(x) that 
satisfies it. A solution that is written in the form ‘y = function of z’ is 
an explicit solution; otherwise it is an implicit solution, i.e. an 
equation of the form F(x,y) = 0 for some function F. 


The general solution of a differential equation is the collection of all 
of its solutions. It is usually possible to give the general solution of a 
first-order differential equation as a formula containing one arbitrary 
constant. A particular solution of a differential equation is a single 
solution containing no arbitrary constant. 


An initial condition for a first-order differential equation 

dy/dx = f(x,y) is an assignment of a value yo that the dependent 
variable y must take when the independent variable x takes some 
given value xo. An initial condition may be specified in the form 
‘y = yo when x = 29’ or ‘y(x9) = yo’; Xo and yo are referred to as 
initial values. 


An initial-value problem is to find the particular solution of a 
differential equation that satisfies a given initial condition. Given the 
general solution of a differential equation involving an arbitrary 
constant C, we can determine the solution satisfying an initial 
condition by substituting the initial values into the general solution; 
this gives an equation from which the required value of C can, in 
principle, be found. 
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5. 


Some differential equations have an analytic solution, i.e. an explicit 
general solution derived using calculus. To decide whether an equation 
of the form dy/dx = f(x,y) may be solved analytically by one of the 
methods described in the module, proceed as follows. 

(a) If f(x,y) is independent of y, so that f(x,y) = f(x), then the 
equation may be solved by direct integration: its general 
solution is y = F(x) +C, where F(z) is an integral of f(x) and C 
is an arbitrary constant. 

(b) If f(x,y) has the form f(x,y) = g(x) h(y), then use the method of 
separation of variables described below. 

(c) If f(x,y) has the form f(x,y) = h(x) — g(x) y (so that the 
equation is of the form dy/dx + g(x) y = h(x)), then the equation 
is linear, and may be solved by the integrating factor method 
described below. 


To solve the differential equation dy/dx = g(x) h(y), where h(y) 4 0, 
using the method of separation of variables, proceed as follows. 


(a) Divide both sides of the equation by h(y), and integrate with 
respect to xz, to obtain 


lien [ae 


(b) If possible, perform the two integrations, obtaining an implicit 
form of the general solution, which should include one arbitrary 
constant. 


(c) If possible, rearrange the formula found in Step (b) to give y in 
terms of x; this is the explicit general solution of the differential 
equation. In addition, there may be supplementary solutions if y 
is a constant satisfying h(y) = 0. 


A first-order differential equation is linear if it can be expressed in the 
form 


W 5 o(x)y = h(x). 


It is homogeneous if h(x) = 0 for all 2, inhomogeneous otherwise. 
For all the cases that you will meet in this module, the initial-value 
problem 
dy 
ae 1 It) Y= W(x), y(z0) = yo, 
ct 
has a unique solution. 


Unit 3 summary 


Linear first-order differential equations can be solved by the 
integrating factor method. An integrating factor for the 
equation dy/dx + g(x) y = h(2) is a function p(x) with the property 
that, after multiplication by p(x), the left-hand side of the equation 
becomes an exact derivative: 


dy d 
n(e) ($4 + a(2)v) = (2) 0). 
Multiplication by an integrating factor therefore makes it possible to 


solve the equation by direct integration. 


The integrating factor for the equation dy/dx + g(x) y = h(x) is 


ple) = exp ( | a(e) ar) 


In the constant-coefficient case dy/dx + Ay = h(x), where A is a 
constant, p(a) = exp(Az). 


The steps of the integrating factor method are as follows. 


(a) Determine the integrating factor p(x). 


(b) Rewrite the equation as (Plz) y) = p(x) h(z). 


(c) Integrate to obtain p(x) y= [re h(x) dx + C, where C is a 
constant. 


(d) Divide through by p(x), to obtain the general solution in explicit 
form. 


Unit 3. Second-order differential equations 


iL, 


A second-order differential equation is a differential equation that 
contains the second derivative d?y/dx? of the dependent variable y 
with respect to the independent variable x, but no higher derivative. 
Thus in addition to d?y/dx?, any or all of the following may occur in 
the equation: x, y, dy/dx. The general solution of such an equation 
normally involves two arbitrary constants. 

A linear second-order differential equation is one that can be written 
in the form 


2, 
a(x) £4 + o(2) © + e(0)y = F(a), 


where a(x), b(x), c(x) and f(x) are given functions. If f(x) is 
identically zero, then the equation is homogeneous; otherwise it is 
inhomogeneous. 
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A linear constant-coefficient second-order differential equation is 
one in which the functions a(x), b(a) and c(x) are constant (with 
a#0), ie. 


a—> +b—+cy= f(z), 
Aim 


and f(x) is a given function. 


The methods of solution described in the unit make use, explicitly or 
implicitly, of the following principle of superposition. 


If yi(x) is a solution of ay” + by’ + cy = fi(x), and yo(z) is a solution 
of ay” + by’ + cy = fo(x), then for any constants k; and ko, 


y(x) = ky yi (x) + ke yo(x) 


is a solution of the equation 


dy | dy 
Te Hos + cy = ky fix) + ke fo(z). 
The general solution of the homogeneous equation 
a di 
a4 4 pF 4 cy=0 


A dx 

is obtained by following the procedure below. 
(a) Solve for \ the auxiliary equation 

a” +bA+e=0. 
(b) (i) If the auxiliary equation has distinct real roots A, and Ao, 

then the general solution of the differential equation is 

y(a) = Ce™* + De®*. 
(ii) If the auxiliary equation has equal real roots A; = Ag 


(i.e. aA? + b\ +c is a perfect square), then the general 
solution of the differential equation is 


y(x) = (C + Da)e™*. 


(iii) If the auxiliary equation has complex conjugate roots 
Ay = a+ if and A2 = a — if, then the general solution of the 
differential equation is 


y(«) = e*(C cos Bx + Dsin Br). 
In each case, C and D are arbitrary constants. 


For the oscillation represented by x(t) = Asin(wt + @), the constant 
A> 0 is called the amplitude, the constant w is called the angular 
frequency, and the constant @ is called the phase constant. The 
period of the oscillation is 27/w. 


Unit 3 summary 


The general solution of the inhomogeneous equation 


is given by y = ye + Yp, where: 


e Yc, the complementary function, is the general solution of the 
associated homogeneous equation 
@y dy 
a—> +b—+cy=0 
dx? " dx avid 
© ¥Yp, a particular integral, is any particular solution of the 
original inhomogeneous equation. 


This leads to the following procedure for finding the general solution 
of an inhomogeneous equation. 


(a) Find the complementary function, by solving the auxiliary 
equation of the associated homogeneous equation. 


(b) Find a particular integral, as described below. 
(c) Add the particular integral to the complementary function. 


For certain inhomogeneous equations, particular solutions can be 
found by the method of undetermined coefficients. This method 
works when the function f(:) is a polynomial, an exponential function 
or a sinusoidal function, i.e. a linear combination of a sine and a cosine. 


The method involves using a trial solution y(x) of a form similar to 
that of f(a), but with coefficients that initially are undetermined. The 
coefficients are determined by substituting the trial solution into the 
differential equation and choosing the coefficients so that the equation 
is satisfied. 


Suitable trial solutions y(a) for specified right-hand-side (or target) 
functions f(a) are shown in the following table. 


Target function f(r) Trial solution y(2) 
Mnx” + My—12—! +++ Pn&” + Pra") +-+- 
+ myx + ™9 + Pix + Po 
mek peke 
moos ka + nsin kr pooskx + qsinkr 
The coefficients in the left-hand column, m,k,mo,m1,..., are given 
constants. Those in the right-hand column, p, po. pi. --., are constants 


to be determined, by differentiating the expression for y(x) twice, 
substituting into the left-hand side of the differential equation, and 
equating coefficients of corresponding terms on the left- and 
right-hand sides. 

Note that even if some of the given coefficients in f(x) are zero — for 
example f(x) = 3x or f(x) = 2cosz ~ it is still necessary in general to 
use the full expression from the right-hand column of the table. 
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10. 


An initial-value problem for a second-order differential equation is 
a problem in which one has to find the particular solution y = y(x) of 
a given equation such that y and its derivative y/ take specified values 
yo and 29, respectively, when the independent variable x takes the 
value zp. The number: ; yo and 20 are called initial values. The 
relationships between initial values are called initial conditions. 
These may be specified either as ‘y = yo and y’ = zo when x = x0’, or 
as ‘y(x0) = yo, y'(0) = 20"- 

An initial-value problem has a unique solution, which can be obtained 
by finding the general solution of the differential equation and 
substituting in the initial values to determine the two arbitrary 
constants that it contains. It is necessary first to differentiate the 
general solution to apply the initial condition y/(29) = zo. 


In a boundary-value problem, a condition is placed on the value of 
either y or its derivative, or some combination of the two, at each of 
two different values of x. The conditions are referred to as boundary 
conditions, and values of x and y in these conditions are boundary 
values. Such a problem may have a unique solution, or no solution, or 
an infinite number of solutions. 


Unit 4 Vectors and matrices 


A vector is a mathematical object consisting of a non-negative real 
number called its magnitude, and a direction. 


In printed text, vectors appear in bold; in handwritten work they are 
written underlined. 


e Two vectors are equal if and only if they have the same 
magnitude and direction. 


e The magnitude of a vector v is denoted by |v| or sometimes by v. 


e The zero vector 0 has magnitude zero; no direction is defined for 
it. 


e The displacement from a point P to a point Q is represented by 
the displacement vector PQ. 


In discussions involving vectors, the word scalar is used to denote a 
real number (positive, negative or zero). 


Scaling a vector or scalar multiplication of a vector is the process 
of multiplying a vector by a scalar. 


For a vector v and (non-zero) scalar m, the scalar multiple mv is 
the vector whose magnitude is |m| |v|, and whose direction is 


e in the same direction as v if m > 0 


e in the opposite direction to v if m < 0. 


Note that —v = (—1)v is the vector with the same magnitude as v 
but pointing in the opposite direction (such vectors are called 
antiparallel). 


Also note that Ov = m0 = 0 for any vector v and scalar m. 

A unit vector is a vector whose magnitude is 1. For v # 0, the 
vector V = (1/|v|)v is a unit vector in the same direction as v. Unit 
vectors along the Cartesian axes are called Cartesian unit vectors, 
and are denoted by i, j and (in three dimensions) k. 

Vector addition is defined geometrically by the triangle rule (see 
the diagram in the margin). 

The vector a+ b is called the sum or resultant of a and b. Vector 
subtraction is defined by a— b = a+ (—b). 

Note that a— a= 0 and 0+ a =a for any vector a. 


A three-dimensional Cartesian coordinate system consists of three 
mutually perpendicular axes, usually labelled x, y and z, that meet at 
a point called the origin. 


A right-handed coordinate system is one satisfying the 
right-hand rule: 


e Point the straightened fingers of your right hand in the direction 
of the positive x-axis, and rotate your wrist until you find that 
you can bend your fingers in the direction of the positive y-axis. 

e Extend the thumb of your right hand. This is the direction of the 
positive z-axis for a right-handed coordinate system. 

Given a Cartesian coordinate system, any vector a can be written 

uniquely in component form as 


a=a,;i+a,j+azk, or equivalently a= (az,ay,az), 
where i, j and k are unit vectors in the directions of the positive z-, y- 
and z-axes, respectively. The numbers az, a, and a, are called the 
(Cartesian) components of a in the directions of i, j and k. The 


process of finding the Cartesian components of a vector is called 
resolving a vector into its components. 


Let a = ari+ ayj + a-k and b = bi + byj + bzk. Then: 
© jaj=/a+az+a2 

e atb= (az +bz)it (ay + by)j + (a: +b.)k 

© ma=(ma,;)i+ (may)j + (maz)k. 


Unit 4 summary 


a+b 


a 


Triangle rule 


Here and below, we give the 
formulas for the 
three-dimensional case; the 
corresponding two-dimensional 
versions are obtained by 
omitting the third component. 
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10. 


The position vector of a point A is the displacement vector OA, 
where O is the origin. So if A has coordinates (x, y, z), then its 


position vector is oA =2i+ yj+ 2k. 
If P and Q are two points with coordinates (xp, yp, zp) and 
(x@, YQ: 2Q); respectively, then the displacement vector POi is 


P@ = 0G - OP = (tq — xp)it (yg — yp)j + (q — 2p)k. 
The vector equation of the straight line joining points P and Q 
whose position vectors are p and q, respectively, is 
r(t)=p+t(q—p)=(1-‘)pt+ta. 
If t is allowed to vary over all real numbers, then r(t) traces out the 


entire straight line through the points P and Q. If 0 <¢ < 1, then r(t) 
traces out the line segment from P to Q. 


The scalar product (or dot product) of two vectors a and b is the 
scalar defined by 


= |al |b] cos, 


where @ (for 0 < @ < z) is the angle between the directions of a and b. 
In particular, non-zero vectors a and b are perpendicular if and only if 
a-b=0. Also, a-a= |al?. 


e = The scalar products of the Cartesian unit vectors are 
i-i=j-j=k-k=1, i-j=j-k=k-i=0 
e In component form, if a = a,i+ a,j +azk and 
b = bri + byj + bzk, then 
a-b= arb; + ayby + azbz. 
The formula |a|? = a- a= a2 + a; +a? is a particular case. 
e The angle @ between two (non-zero) vectors a and b is given by 
Arby + Ayby + azb- 


5 aie /a2 + a3 + a2 [v2 +03 +02 


The component of a vector a in the direction of a unit vector 0 is 


cos 9 = 


ay, =a-t= |alcosé, 
where @ is the angle between a and U. 
The Cartesian components of a are a- i, a-j and a-k; in other words, 


a=(a-i)i+(a-j)j+(a-k)k. 


11. 
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The vector product (or cross product) of two vectors a and b is 
the vector defined by 


aX b= (|a||b/sin@) n, 
where @ (for 0 < @ < z) is the angle between the directions of a and b, 
and ni is a unit vector at right angles to both a and b, whose sense is 
given by the right-hand rule for vector products: this is like the 
right-hand rule for a coordinate system given above, but with the z-, 
y-, z-axes replaced by a, b, ni, respectively. 
The vector product has the following algebraic properties: 
e axb=-(bxa). 
e ax (b+c)=(axb)+(axc) and 

(a+b) xc=(axc)+(bxc). 
e For any scalar A, (Aa) X b = X(a X b) = a x (Ab). 
e ax b=0 if and only if one of a and b is a scalar multiple of the 


other — in other words, one of a and b is zero, or they are parallel 
or antiparallel. (In particular, a x a= 0.) 


e In general, a x (b X c) # (aX b) Xe. 
e The vector products of the Cartesian unit vectors are 

ixi=jxj—kxk=0, 

ixj=—k, jxk=i, kxi=j. 
e Incomponent form, the vector product of the vectors 

a= zi + a,j + ak and b = bi + byj + b-k is 

a X b= (ayb, — azby)i+ (azbr — arbz)j + (arby — aybr) k. 
The area of a parallelogram whose sides are defined by vectors a 
and b is |a x bj. 
The area of a triangle two of whose sides are a and b is sla x bl. 


The volume of a parallelepiped whose sides through one vertex are 
a, b and ¢, is |(a x b) - e}. 

A matrix is a rectangular array of elements (usually numbers) 
arranged in rows and columns. If a matrix has m rows and n columns, 
then it is said to be an m x n matrix (read as ‘m by n matrix’), or to 
be of order m x n. 


e Anm xn matrix A that has a;; as its element in the ith row and 
the jth column is written as 


a1 4412 din 

a1 a22 @2n 
A= 3 3 e 

GQm1i Am2 Amn. 


or A = [aj;] for short. 


Unit 4 summary 


There is an easy way to 
remember this formula in terms 
of determinants: see items 23 
and 27. 
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The second form here is simply 
an alternative to the first, used 
to save space. (The superscript 
T stands for ‘transpose’, which 
is discussed in item 19.) 


36 


14. 


15. 


16. 


e =A matrix with one column is often referred to as a column 
vector; a matrix with one row is a row vector. For example, a 
vector such as a = aji+ a2j + a3k may be written as the column 


vector 
ay 
= = T 
a= |a2 or a=[a, a2 a3)’. 
a3 


e A square matrix is one that has the same number of rows as 
columns. 


e The m xn zero matrix 0 is the m x n matrix all of whose 
elements are zero. 


e Two matrices A = [a,j] and B = [b;;] are equal if they both have 
order m x n, and aj; = 6;; for all i= 1,2,...,m and 


If A = [a;;] and B = [b;;| have the same order, then their sum is a 
matrix of the same order, given by A + B = [aj + jj]. Thus the 
elements of A + B are obtained by adding together the elements of A 
and B in the same positions. For example, for column vectors we have 


ay by ay +b) 
ag} + |b2} = |a2 + be 
a3 bs az + b3 


Similarly for subtraction. 
The scalar multiple of a matrix A = [a;;] by a number 4) is given by 
AA = [Aajj]. 


Thus \A is obtained by multiplying each of the elements of A by 4. 
For example, for column vectors we have 


ay Aa, 
A Jag] = |Aa2g 
a3 az 


Matrix addition and scalar multiplication obey obvious rules of 
algebra. For any matrices A, B and C of the same order, and any 
scalar A, the algebra is: 


e commutative, A+B=B+A 
e associative, (A +B)+C=A+(B+C) 
e distributive, \(A + B) = AA + AB. 


17. 


18. 


19. 


20. 


Unit 4 summary 


The product of an m x p matrix A = [a;;| and a p x n matrix 
B = [bij] is the m x n matrix C = AB, where C = [ej] is formed 
using the ith row of A and the jth column of B, to give 


P 
Cig = Ye ainbsj = ajibij + ajgba; +--+ + Gipbpj- 
k=1 
We imagine placing the ith row of A on the jth column of B, 
multiplying the pairs of elements together, and then adding the 
products. 


The product of two matrices can be formed only if the number of 


columns of the first matrix is the same as the number of rows of the 


second. For example, the product of the matrix A = F ‘| and the 


d 


column vector r= [x y]" is given by 


_ fa b] fa] _ fax + by 

on ke il (F i eae : 
In general, matrix multiplication is not commutative, so AB may not 
be equal to BA. 
For any matrices A, B and C, of the appropriate sizes so that all 
products can be formed, matrix multiplication is: 
e associative, (AB)C = A(BC) 
e = distributive over addition, A(B + C) = AB+ AC. 


The powers of a square matrix are defined in the obvious way: 
A? = AA, A® = AAA, and so on. 


The transpose A? of a matrix A is the matrix obtained by 
interchanging the rows and columns of A. If we denote the matrix A 
by [ajj] and A? by az, then al, = aj. If A is an m x n matrix, then 
A? is an n x m matrix. 
e For any matrices A and B of the same order, (A?)7 = A and 

(A+B)? =A7+B7. 
e For any two matrices that can be multiplied, 

(AB)? = B’A’; 

note the change in order. 
The identity matrix is a square matrix I whose elements are all zero 
except the diagonal elements, which equal 1. For example, the 2 x 2 
10 
0 1) 


Assuming that all the matrix products make sense, the identity matrix 
satisfies AI= A and IA = A. 


identity matrix is I= 
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Although we are free to use any 
row or column when applying 
Laplace’s rule, it is usual to use 
the first row. 
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21. 


22. 


23. 


24. 


The inverse of a square matrix A is the matrix, usually denoted A7!, 
of the same order such that AA“! = ATA =L. 


Only a square matrix can have an inverse, but many square matrices 

do not have inverses: those that do not are called non-invertible (or 

singular), while those that do are called invertible (or 

non-singular). 

e I=L 

e If A and B are invertible matrices of the same size, then AB is 
invertible and 


(AB)! =BA7; 
note the change in order. 


e A2~x 2 matrix 
ab 
a-[ed 
is invertible if and only if ad — be 4 0. When this condition holds, 
we can evaluate the inverse using the formula 


pes fd ad =p 
A -a|! He 


IfA= (: | , then its determinant is a scalar, denoted by det A or 


a 


d 
det A # 0, and non-invertible (singular) if det A = 0. 
The determinant of the 3 x 3 matrix 


ih given by det A = ad — bc. Thus A is invertible (non-singular) if 


a, a2 a3 
A= bh bo bg 
4 C2 C3 
is 
bob by 6: by b 
detA=a;|} “3}—a2|"? °3) 443/71 
ce 6&3 cy 63 ce 


= aj(b2c3 — bgc2) — a2(bic3 — bgc1) + a3(bic2 — bec1). (3) 


The determinant of an n x n matrix for n > 2 can be evaluated using 
Laplace’s rule. 


Given an n x n matrix A = [a;;], its determinant det A may be 
expanded in terms of the elements in row i and their cofactors Cjj as 


det A = aj1Cjy + aigCj2 + --- + @inCin, 


where Cjj = (-1)**47 Mj, and M;; is the minor obtained by deleting 
row i and column j of the original determinant and forming the 
determinant of what remains. 


This rule can be used to derive equation (3). 


25. 


26. 


27. 


28. 


Unit 4 summary 


All determinants have the following properties. 


e Interchanging any two rows or any two columns of A changes the 
sign of det A. 

e = det(A?) = det A. 

e = =©Multiplying any row or any column of A by a scalar k multiplies 
det A by k. 

e = For any scalar k and any n x n matrix A, det(kA) = k" det A. 


e Adding a multiple of one row of A to another row does not 
change det A. Likewise for columns. 


e Ifa matrix A is invertible, then det(A) 4 0 and 
det(A~!) = 1/det(A). 

e For any two matrices that can be multiplied, 
det(AB) = det(A) det(B). 

To find the inverse of an n x n matrix A = [a;;], do the following. 

e Evaluate det A and confirm that det A 4 0. (If det A = 0, then 
no inverse exists.) 


e Evaluate the cofactor Cj; of each element a;;, using the relation 
C,; = (-1)*9 Mj;, where M;; is the minor obtained by deleting 
row i and column j of the original determinant and forming the 
determinant of what remains. 


e Form the square matrix C = [C;j] composed of the cofactors C;;. 
e Take the transpose of C to obtain the matrix C?. 
e Scale the matrix C7 by 1/det A to obtain the inverse of A: 
1 
-1_ T 
aga ° 


The vector product b x c of vectors b = [bz by b.]? and 
c=[c, cy ¢z]7 can be expressed as a 3 x 3 determinant: 


ijk 
bxc=|be by bz}. 
Cy Cy Cz 


lea=[. ag «7, b= & 6,J"ande=—|a. q «/”. 
Then the scalar triple product a- (b x c) is given by 


@, Gy a: 
a-(bxc)=|br by bz}. 
Cr Cy Cz 


e The volume of the parallelepiped with sides defined by the vectors 
a, b and ¢ is given by the modulus of the scalar triple product, 
and therefore by the modulus of the determinant above. 


e Because interchanging rows of a matrix just changes the sign of its 
determinant, the scalar triple product obeys the cyclic identity 


a-(bxXc)=b-(c x a)=c-(axb). 
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29. Any linear transformation of the plane can be represented by a 


2x 2 matrix A, called the transformation matrix. Its effect on the 
position vector r is to transform it into the matrix product Ar. So if 


A= (: | andr=[x yj", then 


_ fa bd] ja} _ jax +by 
ao= [: | A i Pees . 
So the matrix A transforms the position vector r= [x y]” to the 
position vector [ax +by cx +dy]’. 
e The Cartesian unit vectorsi=[1 0]? andj=[0 1]? are 
transformed into the columns of the transformation matrix A. 
e = The unit square whose sides are i and j is transformed into a 
parallelogram with area |det A|, the modulus of det A. 
A singular (non-invertible) matrix transforms the unit square into 
a line (or a point if A = 0). 


e Performing two successive transformations, represented by first A 
then B, is equivalent to the effect of a single transformation, 
represented by the matrix product BA (note the order). 


e = The inverse A~! of a transformation matrix A reverses the effect 
of the transformation. 


e ©The dilation matrix 
ck 0 
D(x, A) = i | 
is a linear transformation that rescales the plane by « in the 
a-direction and A in the y-direction. 


The inverse of a dilation matrix is another dilation matrix: 
D-'(,A) = D(1/«,1/A). This exists only if « 40 and A 40. 
e The rotation matrix 
R(a) = ies a —sin a 
sina cos a 


is a linear transformation that rotates the plane by an angle a 
about the origin, in the anticlockwise direction. 


The inverse of the rotation matrix R(a) is R~!(a) = R(—a). 


Unit 5 summary 


Unit 5 Linear algebra 


Ls 


An equation involving n variables 71, 2r2,...,2, is said to be linear in 
each of those variables if it can be written in the form 


121 + agx2 + 4343 +--+ + Gntn = b, 
where a, @2,...,@,, and 6 are constants. 


Consider a system of linear equations 


4121 + Qy272 +--+ + AinIn = bi, 
9121 + A99%2 + +++ + A2ntp = bo, 


QniZ1 + An2t2 + +++ + Ann&n = bn, 


where there is the same number n of equations as of unknowns 
@1,%,...,2,. Such a system can be written in matrix form as 


Ax=b. 


Here A is the n x n square matrix A = [{a;;], called the coefficient 
matrix, formed from the coefficients a;;. The vectors x and b are the 
column vectors x = (x, x2 ... tal” and b= [b} by ... bn]. 


Usually the elements of A and b are given, while those of x are 
unknown and to be solved for. 


The augmented matrix of the system is the n x (n + 1) matrix that 
has the coefficient matrix for its first n columns and whose final 
column is the column vector b consisting of the right-hand sides of the 
equations. The augmented matrix is usually written Alb, with a 
vertical bar separating the coefficient matrix A from the 
right-hand-side vector b. 

The leading or main diagonal of a matrix A = {a;;| is the collection 
of elements a;;. For a square matrix it is the diagonal that runs from 
the top left corner to the bottom right corner. 


An upper triangular matrix is a square matrix in which each entry 
below the leading diagonal is 0. A lower triangular matrix is a 
square matrix in which each entry above the leading diagonal is 0. 

A matrix that is upper triangular, lower triangular or both 

(ie. diagonal) is sometimes referred to simply as a triangular 
matrix. 

In principle, any system of equations Ax = b, for A non-singular, can 
be solved by calculating the inverse of A using the methods of Unit 4, 
then setting x = A~'b. However, this method is increasingly 
inefficient for systems of three equations or more. A more efficient 
method is to use Gaussian elimination. 
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Performing a row operation 
really involves three sub-steps: 


42 


write down the plan 
implement the plan 
relabel the rows. 


10. 


Gaussian elimination is an efficient and systematic method of 
obtaining the solution x of a system of linear equations Ax = b. 


To solve a system of n linear equations in n unknowns, with coefficient 

matrix A and right-hand-side vector b, by Gaussian elimination, carry 

out the following steps (if possible). 

(a) Formulation: Write down the augmented matrix Alb, denoting 
its rows by Rj,...,Rn- 


(b) Elimination: Adapt the following row operations as necessary. 


(i) Subtract a multiple of R; from Ro, to reduce to zero the first 
element in the first column below the leading diagonal. 

(ii) Similarly, subtract a multiple of R; from R3,...,Ry to 
reduce to zero all the other elements in the first column below 
the leading diagonal. 


(iii) In the new matrix obtained, subtract multiples of Rz from 
R3,...,R,, to reduce to zero all the elements in the second 
column below the leading diagonal. 


(iv) Continue this process until Alb is reduced to U|c, where U is 
an upper triangular matrix. 


(c) Solution: Solve the system of equations with coefficient matrix U 
and right-hand-side vector c by back substitution — that is, 
solve the nth row for x; use this solution to solve the (n — 1)th 
row for z,—1: continue in this way, working up the rows, at each 
stage substituting values already known. 


Sometimes during the Gaussian elimination procedure, an unexpected 
zero can turn up on the leading diagonal element of a row. This can 
make it impossible to use the row for the elimination procedure. To 
remedy this, interchange the row with one below, then continue as 
normal. 


The system of linear equations represented by the matrix equation 
Ax = bis said to be singular when the coefficient matrix A is 
singular (ie. when det A = 0). Such a system does not possess a 
unique solution but, depending on details, may have no solution or an 
infinity of solutions. When the system has no solution, it is sometimes 
said to be inconsistent. 


If Gaussian elimination produces an upper triangular coefficient 
matrix U in which the final element on the leading diagonal is zero, 
then the back substitution process breaks down and the system is 
singular. In this case, if the final element in the final column is 
non-zero, then the system has no solutions, and if it is zero, then there 
may be an infinite number of solutions. 


The equations have a unique solution if and only if none of the 
elements on the leading diagonal of U is zero. 

The n vectors vj, V2,...,V, are linearly dependent if there are 
numbers 01, Q2,...,Q@m, not all zero, such that 


Q1V1 + Q2V2 + +++ + QnVn = 0. 


11. 


12. 


13. 


Unit 5 summary 


In this case we can express one of the vectors as a linear combination 
of the others. 


The vectors are linearly independent if the only solution of the 
equation 


Q1V1 + Q2V2 + +++ +QnVn = 0 
is @ —@5 = 037= ... = en =O 


Two (non-zero) vectors are linearly independent if they are not 
collinear, i.e. parallel or antiparallel. Three (non-zero) vectors are 
linearly independent if they are not coplanar, i.e. if they do not all lie 
in the same plane. 


The space to which vectors belong is often called a vector space. 
The dimension of a vector space is equal to the maximum number of 
linearly independent vectors that it allows. 


In an n-dimensional vector space, if you have n linearly independent 
vectors V1,V2,---,Wn, then you can express any other vector as a 
linear combination of them: 


V=Civ1 + Cove + ++-+EnVn- 


The set of linearly independent vectors vj, V2,.--, v, is called a basis 
for the n-dimensional vector space. 


An eigenvector of a square matrix A is a non-zero vector v such that 
Av = Av for some scalar 4. The number 4 is the eigenvalue 
corresponding to the eigenvector v. Any non-zero scalar multiple of an 
eigenvector is also an eigenvector with the same eigenvalue. 


e Eigenvalues are distinct if they have different values; otherwise, 
they are said to be repeated. Eigenvectors are said to be distinct 
if they are linearly independent. 

e =©Ifa matrix has n distinct eigenvalues, then the corresponding 
n eigenvectors are all linearly independent. 


e  Ifsome of the eigenvalues of a matrix are the same, then the 
corresponding eigenvectors may or may not be linearly 
independent; further investigation is required. 


If an n x n matrix has n linearly independent eigenvectors 
Wy Vajenng Vn; then these can be used as a basis for n-dimensional 
vectors. 


So any n-dimensional vector v can be written as 
V=C1Vi + CoV2 + °+-+CnVn 


for some scalars ¢),€2,.-., c,. This is called the eigenvector 
expansion of v. 


The eigenvector expansion can be used to show that for (almost) any 
vector v and for large k, A*v is proportional to the eigenvector of A 
that has the eigenvalue with the largest modulus. 
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14. 


15. 


16. 


Af. 


Ann xn matrix A has characteristic equation 
det(A — AI) = 0. 


where I is the n x n identity matrix. The left-hand side of this 
equation is a polynomial of degree n in 4. Its roots are the eigenvalues 
of A. 


In principle, the eigenvalues and eigenvectors of a square matrix A can 
be found as follows. Find the roots of the characteristic equation of A 
to obtain the eigenvalues. For each eigenvalue A, solve the 
eigenvector equations (A — AI)v = 0 to obtain the corresponding 
eigenvector. 


If the eigenvalue is not repeated, then this will determine the 
eigenvector up to a scalar multiple. It is usually sufficient to choose 
some value for the scalar multiple to determine a suitable 
representative eigenvector. 


ab 


To find the eigenvalues of the 2 x 2 matrix A = (‘ }: do the 


d 
following. 


(a) Write down the characteristic equation det(A — AI) = 0. 
(b) Expand this as 


a—X b 
ce d-xX 


| = 2? = (a+ a+ (ad ~ be) =0. 

(c) Solve this quadratic equation to find the two values of A, which 
are the required eigenvalues. 

For the above 2 x 2 matrix A, the eigenvalues are 
A=4(at+d4 (a— dP + 4c). 

To find the eigenvector corresponding to the eigenvalue A, do the 

following. 

(a) Write down the eigenvector equations 

(a—Ajz+ by =0, 
cx + (d—A)y =0. 

(b) This pair of equations typically reduces to a single equation that 
is readily solved for x and y. The eigenvector is given by 
v=[zx y\", with x and y replaced by their solved values. Any 
non-zero scalar multiple is also an eigenvector. 


The trace of a square matrix is the sum of the elements on its leading 
diagonal. The trace of A is denoted by tr A. 


The characteristic equation of a 2 x 2 matrix can be written as 


? —trAA+det A =0. 


18. 


19. 


Unit 5 summary 


If the elements of the above 2 x 2 matrix A are real, then the value of 
the discriminant D = (a — d)? + 4be characterises its eigenvalues: 
e If D>0, then there are two distinct real eigenvalues. 


e §=6If D<0, then the two eigenvalues are complex conjugates of each 
other. 


e If D=0, then the two eigenvalues are identical. 


There are a number of general rules that apply to the eigenvalues and 
eigenvectors of an n x n matrix A. To understand these rules, we need 
some definitions about matrices: 


e Areal matrix is one whose elements are all real. 

e Asymmetric matrix is one that is equal to its own transpose: 
A=AT, 

Eigenvalue rules 

e = The product of the eigenvalues of A is det A. 

e The sum of the eigenvalues of A is tr A. 


e Any complex eigenvalues of a real matrix occur in complex _ 
conjugate pairs, i.e. if \ is a complex eigenvalue, then so is A. 

e «©The eigenvalues of a triangular matrix are the diagonal entries. 

e The eigenvalues of a real symmetric matrix are real. 

e A matrix is non-invertible if and only if at least one of its 
eigenvalues is 0. 

Eigenvector rules 

e Ifa matrix has n distinct eigenvalues, then the corresponding 
n eigenvectors are all linearly independent. 

e For a real matrix, the components of an eigenvector corresponding 
to a real eigenvalue can be chosen to be real. 

e For a real matrix, if \ is a complex eigenvalue with corresponding 
eigenvector v, then V, the vector whose components are the 
complex conjugates of those of v, is an eigenvector corresponding 
to the eigenvalue A. 

e For any real symmetric matrix A, the eigenvalues are real, and 
the eigenvectors vj; may be chosen to be real and orthogonal to 
each other, i.e. vj + vj = 0 for i # j. Using matrices, the inner 


product v; - vj is calculated 
i 2 
using V; Vj- 
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Note that there are exceptions to 
this rule, but we do not consider 
such cases in this module. 
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Unit 6 Systems of linear differential equations 


1. 


2. 


A system of linear constant-coefficient first-order differential 
equations for unknowns 2}, £2,...,2;, takes the form 


Ey = A112] + ay2r2 +--+ + intr + hi(t), 
£2 = A121 + a2Q%2 +--+ + A2nFp + halt), 


Fn = Ani + An2%2 +++ + nn Tn + hy(t), 
where there is the same number n of equations as of unknowns. 


The independent variable t often represents time. The dependent 
variables are x(t), 22(t),..-, x,,(t); but for systems with n = 2 or 
n = 3, x, y and z are often used instead of x1, x2 and x3. The 
coefficients aj; are constants. 


The equations may be expressed in matrix form as 


x=Ax-+h, 
where 
x=[e1 a ... onl’ and x=[%) ao ... an)’; 


here A = [{a;;| is a constant matrix called the matrix of coefficients, 
but the vector h= [hy ho ... hal” may depend on t. 


The general solution of a first-order system of n equations can be 
expressed as a solution containing n arbitrary constants. 


A particular solution is a solution containing no arbitrary constants 
and satisfying given conditions. 

Particular solutions of systems of n first-order differential equations 
are usually obtained by demanding that the general solution satisfies 
n initial conditions. These conditions fix the n arbitrary constants. 

A system with h = 0, i.e. of the form x = Ax, is said to be 
homogeneous. An inhomogeneous system is one where h 4 0. 


To find the general solution of x = Ax, where A is an n x n matrix, 
proceed as follows. 


Tf the matrix A has only real eigenvalues: 
e Find the eigenvalues Aj, A2,..., A, and corresponding eigenvectors 
Vi, V2,---,;Vn of the matrix A. 


e Write down the general solution in the form 


Aut 


x = Cyvie™! + Covoe** +--+ + Cavne, 


where C},C2,..., C,, are arbitrary constants. 


Unit 6 summary 


If the matrix A has some complex eigenvalues (which must occur in 
complex conjugate pairs \ and A, with corresponding complex 
conjugate eigenvectors v and V), then a further step is necessary: 


e Replace the complex terms ve and vert appearing in the general 
solution with Re(ve™) and Im(ve™). 


The general solution will then be real-valued for real 
C1, Co,...,Cn- 


The above procedure fails if A does not have n linearly independent 
eigenvectors. This case is not covered in this module. 


The complementary function of an inhomogeneous system 

x = Ax +h is the general solution of the corresponding 
homogeneous system x = Ax. A particular integral of the 
inhomogeneous system is any solution of it containing no arbitrary 
constants. 


The general solution of the inhomogeneous system x = Ax + h takes 
the form x, + Xp, where x, is the complementary function and x, is a 
particular integral. 


To find a particular integral x, = [xp Yl? of x = Ax +h, where A 

is a 2 x 2 matrix, try a solution constructed as follows. A similar procedure works for 
= “ systems of more than two 

e When the elements of h are polynomials of degree k, choose xp equations. 


and yp to be polynomials of degree k. 


e When the elements of h are multiples of the same exponential 
function, choose zp and y, to be multiples of this exponential 
function. 


The trial solution will have a number of undetermined coefficients. 
To find their values, substitute the trial solution into the system of 
differential equations and equate coefficients of corresponding terms. 
This will give a number of equations that can be solved for the 
coefficients. 


A system of homogeneous linear constant-coefficient 


second-order differential equations for unknowns 2}, £2,...,2p 
takes the form 


#, = ay2, + ayo%2 +--+ Aintn, 


Hq = ag 21 + A22%2 + +++ + Ann, 


Ep = Gn1@1 + Gn2te +++: +Gnn2n, 
where there is the same number n of equations as of unknowns. 
The equations may be expressed in matrix form as 

* = Ax, 


where A is a constant matrix. 
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8. 
Complex eigenvalues and 
repeated real eigenvalues are not 
discussed, but they can be dealt 
with by generalising here. 

9. 


To solve a system X = Ax, where A is an n x n matrix with n distinct 
real eigenvalues, do the following. 


e Find the eigenvalues Aj, A2,..-, An of A, and a corresponding set 
of eigenvectors V1, V2.---,Vn- 


e Each positive eigenvalue 4, corresponding to an eigenvector v, 
gives rise to two linearly independent solutions 


Vit Vx 


ve and ve 


Each negative eigenvalue \, corresponding to an eigenvector v, 
gives rise to two linearly independent solutions 


vcosV—At and vsinV—At. 


A zero eigenvalue corresponding to an eigenvector v gives rise to 
two linearly independent solutions 


v and vi. 


e = The general solution is then an arbitrary linear combination of the 
2n linearly independent solutions found in the previous step, 
involving 2n arbitrary real constants. 


Systems of second-order differential equations arise in the study of 
oscillating systems. A normal mode of oscillation is one in which 
all of the coordinates of the system oscillate sinusoidally with the same 
angular frequency. 


Unit 7 Functions of several variables 


In this module, unless told 
otherwise, you may assume that 
functions are ‘sufficiently 
smooth’ for their derivatives to 
be defined and well-behaved. 
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A function f(x,y) defines a surface z = f(x,y), which can be 
visualised in a perspective view. If one of the variables is given a fixed 
value (y = a, say), then the corresponding section function f(:,a) 
can be plotted as a normal graph. A contour map for f(x,y) is a 
collection of contour lines in the ry-plane. Along each contour line, 
f(x,y) has a constant value. 


For a function f(2,y,z) of three variables, the contour lines are 
replaced by contour surfaces. 


The partial derivative Of /0x of a function f(x,y) is obtained by 
differentiating f(x,y) with respect to x, treating y as a constant. It is 
also written as f,(x,y) or f,. All the usual rules of differentiation 
apply. More formally, 


oa = frla,y) = Yim LE +89) = fey) 


6x0 ox 
Similar definitions apply to Of /Oy = f,(x,y), and to partial 
derivatives of functions of more than two variables. In a partial 
derivative with respect to a given independent variable, all the other 
independent variables are treated as constants. 


Unit 7 summary 


Functions obtained by partially differentiating a function once are 
called first-order partial derivatives. These can be partially 
differentiated again to obtain higher-order partial derivatives. For 
example, 


Pf of Of of 
Ox -z( 3) = fee and Oy Ox = alas a= i= 


are second-order partial derivatives. 


For a sufficiently smooth function f(x,y), the mixed partial 
derivative theorem states that 


of _ of 

Oxdy Oy Ox’ 

A similar result applies to a function f(21,22....) of n variables. 
Again, the order of differentiation does not matter, so 


or equivalently fry = fyz- 


for all a; and x;. 


The chain rule of partial differentiation takes several different 
forms. 


e = Chain rule for small changes If f = f(x,y), and the independent 
variables « and y change by small amounts dz and dy, then the 
corresponding small change in f is 

or PF sy 
afia 

FS Oy 

e =Chain rule “a differentiation with respect to a parameter If 
f = f(x,y), and the independent variables x = x(t) and y = y(t) 
are functions of a parameter ft, then the rate of change of f with 
respect to t is 


df _ af dx , af dy 


dt Ox dt ay dt’ 


e = Chain rule for a change of variables If f = f(x,y), and the 
independent variables x = z(u,v) and y = y(u,v) are functions of 
two other variables, u and v, then the partial derivative of f with 
respect to u is 


Of _ of On | Of Oy Oy 


Qu Ax Ou | Oy Ou" 


with a similar result for the partial derivative of f with respect 
to v. 


et a 


All these chain rules can be extended to functions of more than two 
variables. For example, if f = f(x,y,z), then the chain rule for small 
changes becomes 


apes as n+ Shy vt SE 0 
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Given a surface z = f(x,y), and a direction in the ry-plane defined by 
a unit vector N = cos4i-+ sin @j, the slope of the surface in the 
direction of the unit vector is 


of 
lope = cos@ 
slope = cos 0 


<pOt x ge 
Oa pant =z fr + iy fy- 


For a function f(x,y), the gradient vector is vector-valued function 
of x and y in the xy-plane given by 


It has the following properties: 


e For a unit vector i in the ry-plane, the slope of the surface 
z= f(x,y) in the direction of the unit vector is 


slope = fi- grad f. 
e grad f points in the direction in the ry-plane that gives maximum 
slope, and its magnitude is equal to the maximum slope. 


e At any given point, grad f is perpendicular to the contour line 
through the point. 


The first-order Taylor polynomial for f(x,y) about (a,b) is 
pi(x, y) = f(a.) + fr(a,b)(%— a) + fy(a,b)(y — 6). 

This matches the values and first-order partial derivatives of f(x,y) 

at (a,b). 

The plane given by 


z=pilz,y) 


is called the tangent plane for f(x,y) at (a,b). For any direction in 
the xy-plane, the slopes of the function f(x,y) and the corresponding 
tangent plane are identical at (a,b). 


The second-order Taylor polynomial for f(x,y) about (a,b) is 


pa(x, y) = pila. y) 
+ 3(fex(a,b)(x — a)” + 2fey(a,)(x — a)(y —b) + fyy(asb)(y — b)”) - 

This matches the values and first- and second-order partial derivatives 
of f(x,y) at (a,b). 
A point (a,b) inside the domain of f(x,y) is a stationary point if 
both f,(a,b) and f,(a,b) are equal to zero. There are three types of 
stationary point. 
e =A local minimum occurs at (a,b) if there is a small region 

around (a,b) within which f(x,y) > f(a,6) at all points 

(x,y) # (a,b). 
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e =A local maximum occurs at (a,b) if there is a small region 
around (a,b) within which f(x,y) < f(a,b) at all points 
(x,y) # (a,b). 

e A saddle point is a stationary point that is neither a local 
minimum nor a local maximum. Through such a point, some 
paths climb to higher function values while others descend to 
lower function values. 


Stationary points are found by taking all the first-order partial 
derivatives of the given function, setting them equal to zero, and 
solving the resulting set of simultaneous equations. For example, the 
stationary points of f(x,y) occur at points (a,y) that satisfy the 
simultaneous equations 
fee.y)=0 and f(x,y) =0. 
The Hessian matrix is the matrix of the second-order partial 
derivatives of f. 
To classify the stationary point at (a,b) of a smooth function f(x,y), 
do the following. 
(a) Find the second-order partial derivatives, and evaluate them at 
the stationary point. 
(b) Construct the Hessian matrix H at the stationary point, and find 
its eigenvalues. 


(c) Apply the following rules: The test is inconclusive if any of 


é yc the ei ues are equal to zero. 
e ‘If all the eigenvalues are positive, then we have a local eae een eae 


minimum. 


e If all the eigenvalues are negative, then we have a local 
maximum. 


e If the eigenvalues have mixed signs, then we have a saddle 


point. 
Alternatively, for a function f(x,y) of two variables, the determinant 
test can be used. In this case the Hessian matrix can be written as The Hessian matrix of a smooth 
Ges [ee b) fry(a, Hy _ z | function is always symmetric. 
Fiye(a,b)  fyy(a,b) iB iC. 
The stationary point is: 
e a local minimum if det H = AC — B? > 0 and A>0 We could equally well use C 


instead of A in the determinant 
test. 


e a local maximum if det H = AC — B? > 0 and A <0 
e asaddle point if det H = AC — B? < 0. 


The determinant test is inconclusive if det H = 0. 
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Unit 8 Multiple integrals 


1. 


The area integral of a function f(x,y) over a rectangular region S in 
the ry-plane, bounded by the lines z = a, rx = b and y=c, y =d, can 
be found by two successive integrations. The inner integral is always 
performed first. For example, we can write 


[senaa- [ : ( [" Heal ar) em 


where the integral over z is carried out first, with y held constant. 
This gives a function of y, which is integrated over y. 


Alternatively, the same area integral can be written as 


[senaa= [ ( ia te.) ay) de. 


In this case, the integral over y is carried out first, treating x as a 
constant; the result is then integrated over x. 


The volume integral of a function f(x,y, z) over a cuboid region R, 
bounded by the planes « = a), x = a2, y= bi, y = bo, z = cy and 

z =», can be found by three successive integrations. The innermost 
integral is evaluated first, and successive integrations move 
progressively outwards. For example, we can write 


[ fewaav = [ im (f° s2,2) 42) au) dx, 


where the integral over z is done first, with x and y held constant. 
The result is integrated over y, with x held constant, and the final 
integral is over x. 


Alternative orderings may be chosen. Inner integrals are always 
completed before integrals placed outside them. The limits of 
integration in a given definite integral can depend only on the variables 
of integration in integrals that lie further outside it (and are done 
after it). The limits of the outermost integral are always constants. 

In an area integral, if the function to be integrated takes the product 
form f(x,y) = g(x) h(y) and the limits of integration are all 
constants, then the area integral can be expressed as a product of 
definite integrals over x and y: 


[senaa= [ ogee [ =~ aly) dy. 


=a =e 
In a volume integral, if the function to be integrated takes the product 
form f(x,y, 2) = u(x) v(y) w(z) and the limits of integration are all 
constants, then the volume integral can be expressed as a product of 
definite integrals over w, y and 2: 


| few. z)dV= "es u(x) dx x i- u(y) dy x i w(z) dz. 


However, such factorisations do not work in general. 


Unit 8 summary 


To find the area integral of a function f(x,y) over any region S in the 

ay-plane, start by choosing which integral to do first — that over x or 

that over y. The following steps assume that the integration over y is 
done first. 

(a) Draw a diagram showing the region of integration S. 

(b) Draw a vertical strip parallel to the y-axis, centred on x, and 
spanning the region (as shown in the diagram). Determine the 
lower limit y = a(x) and the upper limit y = §(x) for this strip. 
These are the limits for the y-integration (the inner integration). 
In general, they are non-constant functions of x. 


(c) Determine the minimum value x = a and the maximum value ¢ w oe 
ax = b for x-values throughout the region. These are the limits for 
the z-integration (the outer integration), and are always 
constants. 


Sketching a region of 
integration when the integral 
over y is done first 

(d) Write down the area integral as 


I few)da = [ mn ( i ine steau) de. 


(e) Evaluate the inner integral over y first, holding 2 constant, and 
substitute in the limits of integration. This gives a function 


=A) 
g(x) = | J f(x,y) dy. 
y 


=a(r) 
(f) Evaluate the remaining definite integral of g(x) over x. 


If you choose to do the z-integral first, then the limits of integration 
must be found using a different sketch, with horizontal strips running 
parallel to the z-axis. The lower and upper limits of the z-integration 
are then functions x = u(y) and x = v(y), and the lower and upper 
limits of the y-integration are constants y = c and y = d. The integral 
is then written in the form 


feoapaae [| Eset) ay. 
| ie a=u(y) 


The integral over x is done first, with y held constant; the result is 
then integrated over y. 


The volume integral of a function f (x,y,z) over any region R can be 
written as 


[semsar= [fOr (fe seats) a) a 


The limits of integration can be found by drawing two diagrams — a 
perspective view of the region R, and a projection of R onto the 
zy-plane. The limits of the 2-integration are found from the 
perspective view, by considering a thin column centred on a typical 
point (2,y). The limits of the y-integration (which are functions of x) 
are found by considering a typical thin slice parallel to the y-axis, 
represented by a thin strip in the projection onto the ry-plane. 
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The constant limits for the outer z-integration are also found from 
this projection. This method can be adapted for different orders of 
integration. 


The area of a region in the xy-plane is equal to the area integral of the 
function f = 1 over the region. Similarly, the volume of a region in 
three-dimensional space is equal to the volume integral of the function 
f =1 over the region. 


Polar coordinates (r, 0) are related to Cartesian coordinates (x, y) 
by the equations 


“=rcosd, y=rsing. 
In polar coordinates an area element has area 
6A =r or dd. 


The area integral of a function f(r,@) over a disc S of radius R 
centred on the origin is 


[seoaa= [ - ( | i f(r.d)r ar) de, 


where we have chosen to integrate first over r, and then over ¢. The 
reverse ordering is also valid. The limits of integration can be adjusted 
for regions occupying parts of a disc. 

Cylindrical coordinates (r,@, 2) are related to Cartesian 
coordinates (x,y,z) by the equations 


@=rcos?, y=rsng, z=2. 


In cylindrical coordinates a volume element has volume 


OV =r or b¢6z. 


aL 


(a) Cylindrical coordinates; (b) a volume element in cylindrical coordinates 


(b) 


Unit 8 summary 


The volume integral of f(r, @,z) over a cylindrical region D, centred on 
the z-axis with radius R, and with flat surfaces at z = 0 and z = h, is 


I f(r,4,2) dV = [. (f ([. f(r,6,2) rar) a) dz. 


Alternative orderings are also valid. The limits of integration can be 
adjusted for regions based on parts of a cylinder. 


9. A region that is axially symmetric around the z-axis has volume 


Van f” (rhac(2) ~rbin(2)) de 


=21 


where Tmin(Z) and rmax(z) are the minimum and maximum radii of the 
region for a given value of z, and z; and z2 are the minimum and 
maximum values of z. For a solid object with no holes, ryin(z) = 0. 


10. Spherical coordinates (r,6,@) are related to Cartesian coordinates 
(x,y,z) by the equations 


xz=rsinOcosé, y=rsinésing, z=rcosé. 
In spherical coordinates a volume element has volume 
5V =r? sin dr 6060. 


ay 


aa 
BD 


R 


(a) (b) 
(a) Spherical coordinates; (b) a volume element in spherical coordinates 


The volume integral of f(r,@,@) over a spherical region of radius R, 
centred on the origin, is given by 


o=2n =n r=R 4 
I= i ( f) ( | f(r,0,6) r° sind dr) a) do. 
o=0 0=0 r=0 


Alternative orderings are also valid. The limits of integration can be 
adjusted for regions based on parts of a sphere. 

11. A coordinate line is a line along which one coordinate varies while 
the other coordinates remain fixed. 
A coordinate system is said to be orthogonal if its coordinate lines, 
for different coordinates, meet at right angles. Cartesian, polar, 
cylindrical and spherical coordinate systems are all orthogonal. 
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12. 


13. 


For any coordinate u, the length of the segment of the u-coordinate 
line between u and u+ du, where du > 0, is expressed as 


length of segment = h, du, 


where h,, is called the scale factor for the u-coordinate; this may be a 
function of the coordinates. 


Tf the coordinates (u,v, w) are related to Cartesian coordinates by 
equations of the form 


r=x(u,v,w), y=y(u,v,w), = 2(u,v,w), 


then the scale factor h, is given by 


n= (2) +2) +) 


with similar formulas for h, and hy. 


Scale factors in some orthogonal coordinate systems are given in the 
table below. 


Cartesian (x, y) hy =1, hy =1 

Cartesian (z,y,z) he=1, hy=1,h,=1 
Polar (r,@) h,=1,hg=r 

Cylindrical (r,¢,z) hy =1,hg =7r, hz =1 
Spherical (r,6,¢) hy =1, ha =r, hg =rsind 


In any orthogonal coordinate system (u,v) with scale factors h, 
and h,, an area element has area 


6A = hyhy du dv. 


In any orthogonal coordinate system (u,v,w) with scale factors hy, hy 
and h,,, a volume element has volume 


OV = hyhyhy du dv bw. 


A surface can be parametrised by two coordinates (u,v). These are 
related to Cartesian coordinates by equations of the form 


v=a2(uv), y=y(u,v), z=2(u,v). 


Using these equations, we define the Jacobian vector 


ijk 
Ox Oy Oz 

J=|du du dul: 
Ox Oy Oz 
Ov dv av 


which is perpendicular to the surface. 


Unit 9 summary 


Then a surface area element has area 
5A = |J| Ou dv, 
where |J| is the magnitude of J. 


For a surface region S with uw; <u < ug and v < v < ve, where the 
limits of u and v are all constants, the surface integral of a function 


f(u,v) is 


[sea=[- in f(u,v) is\du) dv. 


For f = 1, this gives the area of the surface region S. If the surface S 
is part of the surface of the sphere of radius R centred at the origin, 
and (u,v) are the spherical coordinates (0, ¢), then J = R? sine,. If 
S is part of a cylinder of radius R with axis along the z-axis, and 
(u,v) are the cylindrical coordinates (@, z), then J = Re,. 


Unit 9 Differentiating scalar and vector fields 


Ts 


A field is a quantity with definite values at points throughout a region 
of space. Scalar fields and vector fields describe the distribution of 
scalar and vector quantities, respectively. At any given point: 


e ~The value of a scalar field is independent of the orientation of the 
coordinate system. 


e The magnitude and direction of a vector field are independent of 
the orientation of the coordinate system. (This implies that the 
components of a vector field depend on the choice of coordinate 
system.) 


In a coordinate system (u,v, w), the unit vector e, associated with 
the coordinate u is a vector of unit length pointing in the direction in 
which wu increases while the other coordinates v and w are held 
constant. In orthogonal coordinate systems, the unit vectors are 
mutually orthogonal. They may vary from point to point. 


In Cartesian coordinates (x,y,z), the unit vectors are 
ep =i, ey=j, e.=k. 
In polar coordinates (r,@), the unit vectors are 
e, =cosdi+singdj, es = —sindi+cos dj. 
In cylindrical coordinates (r,@,2), the unit vectors are 
e, =cosdit+sindj, eg = —sindi+cos@j, e,=k. 
In spherical coordinates (r,0,@), the unit vectors are 
e, = sin@cos di + sin #sin dj + cos Ok, 
e9 = cos 6 cos di + cos @sin gj — sin #k, 
eg = —sindi+ cos@j. 
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is often used. 
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To convert a vector field F from Cartesian coordinates into another 
orthogonal system (u,v, w), for each component do the following. 
(a) Write down 
Fy = eu+F, 
and expand the scalar product on the right-hand side using 
Cartesian expressions for e,, and F (involving i, j and k). 
(b) The resulting expression generally depends on (x,y,z). Use the 
coordinate transformation equations in the form 
@=a2(u,v,w), y=y(u,v,w), z= 2(u,v,w), 
to obtain an expression for F,, solely in terms of u, v and w. 
Once this has been done for each component, write down 
F = Fyeu + Fy ev + Fw ew. 
Given a scalar field V(x, y, z), the corresponding gradient vector 
field is 
We= qeit pit gk 
This is also called the gradient of V. 
At each point, the gradient has the following properties: 


e Its direction is that in which V increases most rapidly. This 
direction is perpendicular to the contour surfaces of V. 


e Its magnitude is the maximum rate of increase of V with respect 
to distance travelled in three-dimensional space. 


e = The gradient of a scalar field is a vector field, so its magnitude 
and direction are independent of the orientation of the coordinate 
system. 


e Fora small displacement ds = dri + dyj + 62k, the corresponding 
small change in V is 


OV = VV - ds. 
e = The rate of change of V with distance in the direction of the unit 
vector ni is equal to the component of VV in the direction of n: 


rate of change of V = n- VV. 


In polar coordinates, 
ov 10V 
w= -> &. 
or + 5 ag 
In cylindrical coordinates, 


ov 1av ov 


W= a ets ag Ot ae 
In spherical coordinates, 
ov 10V 1 OV 
We ar + ag + rain 06 
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Given a vector field F = F,i+ Fj + F.k, the corresponding 
divergence is 
V-F= oP + a + OF The alternative notation div F is 
Ow oy Oz often used. 
The divergence of a vector field F is a scalar field. It describes the 
extent to which F flows outwards or diverges from each point. 


In polar coordinates, 

1a) , 1 OF 

r Or r do" 

In cylindrical coordinates, 

10(rF,) 10F, , OF: 

Fs 2 Foe | Oz” 

In spherical coordinates, 

1 O(r?F,) 1 Asin Fe) 1 OFs 
rr? Or e rsin@ 00 rsin@ 06° 
Given a vector field F = F,i+ Fj + Fk expressed in a right-handed 
Cartesian coordinate system, the corresponding curl is 


V-F= 


V-F= 


V-F= 


ij k 
ao 0 0 
VxF=|=— = al I The alternative notation curl F 
Ox Oy dz is often used. 
Fi. K F, 


In the determinant, the partial derivative operators in the second row 
act on the components in the third row. Expanding out the 
determinant gives 
OF, OF, OF. OF, OF, OF, 
F= =——4)i- = — j #_—=)k. 
ve (F me) (FE Oz ae Ox Oy 
The curl of a vector field F is another vector field. It describes the 
extent to which F rotates or swirls locally about each point. (The 
sense of the rotation is given by the right-hand rule.) 


In polar coordinates, 


O53) 


In cylindrical coordinates, 


1 
VxF=- 
r 


e, Treg e 


i110 8 98 
VxXE= "5 de Dal: 
Fe rFy Fe 


In spherical coordinates, 
e, reg rsinfdeg 
il (ne) a) 
~ p2sin@|Or 00 do 
F, rFg rsind Fy 
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2. 

3. 
This formula can be extended in 
an obvious way for paths in 
three-dimensional space. 

4. 
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Integrating scalar and vector fields 


A path is a curve with a definite sense of progression from a start 
point to an end point. If the start and end points are identical, then 
the path is closed; otherwise, it is open. 


A path can be represented by a set of parametric equations of the 
form 


r=c(t), y=y(t), 


where x, y and z are Cartesian coordinates of points on the path, and 
the parameter ¢ increases monotonically, from ¢ = ¢; at the start point 
to t = tg at the end point. For a path in the xry-plane, we need only 
the parametric equations for x and y. 


z=2(t) (t1<t<tr), 


A path in three-dimensional space, with parametric equations 
x= x(t), z=2(t) (tr. <t< tr), 
has length 


y=y(t), 


If the path is restricted to the xy-plane, then this simplifies to 


ts dey? dy\? 
a= ['V(Z) +(@) 
The line integral of a scalar function A(x, y) along a path with 
parametric equations x = x(t), y = y(t) for t) < t < tg is given by 


M= [* scat.) (4) +(#) a 


If A(x, y) is the mass per unit length along the path, then the line 
integral gives the total mass along the path. 
To calculate the line integral of a vector field 
F = F,i+ Fyj+ Fk 
along a path C with the parametric representation 
r=a2(t), y=y(t), z=2(t) (h <t< tp), 
do the following. 
(a) Use the parametric representation to find 
ds dx, dy, dz 
ad! at) at 
(b) Express the components of F as functions of the parameter t. 
(c) Find the scalar product 
ds dx dy dz 


Pen at ate 


as a function of t. 
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(d) Evaluate the line integral as a definite integral over ¢: 


te 
[Feas= [P-Sat 
C fi dt 


A vector field F is called a gradient field if it can be expressed in the 

form 
F=-—VU, 

where U is a scalar field, known as the scalar potential field 

associated with F. 

Line integrals of gradient fields have the following properties: 

e Any line integral of a gradient field is path-independent (i.e. it 
does not depend on the detailed path between given start and end 
points). 

e Any line integral of a gradient field around a closed loop is equal 
to zero. 

A vector field F is said to be conservative if, throughout its domain 

of definition, all its line integrals are path-independent. 

The line integral of a conservative field F along a path C from a start 

point A to an end point B may be written as 


| F -ds, 
A+B 


as there is no need to indicate the precise path C. 


For a conservative field F with an associated scalar potential 
field U(r), 


| F - ds = U(r,) — U(rp). 


All gradient fields are conservative, and all conservative fields are 
gradient fields, so the terms conservative field and gradient field are 
synonymous, and can be used interchangeably. 

Given a conservative vector field F, the corresponding scalar potential 
field is defined as follows. We choose a point ro as the position of zero 
scalar potential; by definition, U(r9) = 0. At any point r, the scalar 
potential field is then given by the line integral 


U(r) = -[ F -ds, 
Tor 


where the path may be taken to be a straight line from ro to r, or any 
other choice that simplifies the integration. A useful check is that U(r) 
should satisfy 


F=-VU. 
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In this module all regions are 
simple. 


An oriented area 
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10. 


1. 


The curl test is used to determine whether a given vector field F is 
conservative in a simple region (such as R°). The vector field F is 
conservative if and only if V x F = 0 throughout the region. 


If is a unit vector perpendicular to a tiny planar surface element of 
area 6S, then the oriented area of the element is 


6S = 6Sn. 


Given a vector field F and a planar element with oriented area 6S, the 
flux of the vector field over the element is defined as 


flux = F - dS, 


where the field F is evaluated at the position of the element. This is 
the normal component of the field (i.e. the component in the direction 
of the unit normal) multiplied by the area of the planar element. Flux 
is a scalar quantity that can be positive, negative or zero depending on 
the relative orientations of F and the unit normal n. 


A closed surface is one that divides space into two regions — an 
exterior space and an interior space. An open surface is one that is 
not closed. When a surface is approximated by many small planar 
elements, neighbouring elements are chosen to have unit normals that 
are almost parallel (rather than almost antiparallel). For any closed 
surface, all the unit normals are chosen to point outwards into the 
exterior space, rather than inwards towards the enclosed volume. 
Suppose that a surface S is parametrised by (u,v), for u) <u < ug 
and v, < v < v2, where the minimum and maximum values of the 
parameters are constants. Then the flux of a vector field F over a 
surface S is given by the surface integral 


v=v2 u=u2 
tux= [ F-as = [ (/ F-Jadu) de, 
Ss v=v1 u=uL 


where 
ij k 
Ox Oy Oz 
J=|du du du 
Ox Oy Oz 
Ov dv dv 


is the Jacobian vector in item 13 for Unit 8. To evaluate the surface 
integral, the integrand F - J must be expressed in terms of the 
parameters u and v. 


On the surface of a sphere of radius R, centred on the origin and 
parametrised by (@,@) of spherical coordinates, 


J= R’sinde,, 
where 
e, = sin@cos@i+sin@sin dj + cosdk 


is the radial unit vector of spherical coordinates. 
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13. 


14. 


15. 


16. 


17. 


Unit 10 summary 


The divergence of a vector field F at a given point is related to the 
flux of F over a tiny surface enclosing the point. In the limit where 
the surface area and its enclosed volume shrink to zero, 

V-F-= flux of F over surface 
~ volume enclosed by surface’ 
So the divergence of a vector field at any point can be interpreted as 
the flux per unit volume at that point. 
Given a vector field F and a closed surface S enclosing a volume V, 
the divergence theorem states that 


[F-= [ v-Fav. 
s Vv 


This is also called Gauss’s theorem. 


The unit normal of a planar element and the sense of positive 
progression around the perimeter of the element are related by the 
right-hand grip rule: with the thumb of your right hand pointing in 
the direction of the unit normal of a planar element, the curled fingers 
of your right hand indicate the sense of positive progression around 
the perimeter of the element. 


Given a vector field F and a closed path C, the circulation of F 
around C is given by the line integral 


The right-hand grip rule 


circulation = ¢ F - ds. 
Cc 


The circle on the integral sign is optional: it is used to indicate that 
the path is closed. If C is a closed path around a planar element with 
a given unit normal, then it is understood that C is traversed in the 
positive sense determined by the right-hand grip rule. 

Given a vector field F in the vicinity of a given point, the component 
of V xX F in the direction of the unit vector A can be found by taking 
a planar element with unit normal fi at the point. In the limit where 
the element becomes very small, 


circulation around perimeter of element 
area of element : 


So each component of the curl at a given point can be interpreted as a 
circulation per unit area at that point. 


(V xF)-n= 


If F is a vector field and S is an open surface with perimeter C’,, then 
the curl theorem states that 


fF eds= [ov x¥)-as, 


where it is assumed that C is traversed in the sense determined by the 
right-hand grip rule. This result is also called Stokes’s theorem. 
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Unit 11 Fourier series 


1. 


A function f(t) is periodic if there is some positive number A such 
that f(t +A) = f(t) for all t. The number 4 is a period of the 
function f. If A is a period of a function, then so is nA for any positive 
integer n. The fundamental period of a periodic function is the 
smallest (positive) value for the period. The symbol rt is often used for 
the fundamental period. Unless it is specified otherwise, assume that 
‘period’ means ‘fundamental period’. A fundamental interval for a 
periodic function is any interval whose length is the fundamental 
period. 


The function f(t) = cos(w#) is periodic with fundamental period 27/w. 
The constant w is called the angular frequency. 


A function f(t) is even if 
f(—t) =f(#) for all values of t. 
For any integer k, cos(kt) is even. 
A function f(t) is odd if 
f(—t) =—f(t) for all values of t. 
For any (non-zero) integer k, sin(kt) is odd. 


Even and odd functions combine under addition and multiplication as 
follows. 


e The sum of two even functions is even. 

e The sum of two odd functions is odd. 

e = The product of two even functions is even. 

e = The product of two odd functions is even. 

e The product of an even and an odd function is odd. 


The integrals of even and odd functions can be simplified when the 
range of integration is symmetric. 


e = If f(t) is even, then 
. dt =2 : dt. 
[" sleyae=2 [plea 
e = If f(t) is odd, then 


” Fiat =0. 


—a 


6. 


Unit 11 summary 


A Fourier series is an infinite series of trigonometric functions of the 
form 


Ay + Ay os (2) + Ascos (=) +::- 
i r 
4 
+ Bysin (=) + Basin (=) 42 y 
T T 


where the constants Ag, Aj, A2,..., By), Bo,... are called the Fourier 


coefficients, and 7 > 0. 


A Fourier series defines a periodic function whose fundamental period 
is T. 


To find the Fourier series for a periodic function f(t), do the 
following. 


(a) Find the fundamental period 7. 
(b) Write down the Fourier series 


F(t)= fi S cos (4 


n=1 


UT 7 
)+ 5 B, sin (# ‘), 
T 


n=1 


where Ap and the A, and B,, are coefficients to be determined. 
Simplify the arguments of the sines and cosines where possible. 


(c) Use the following formulas to determine the Fourier coefficients: 


7/2 
Ap = =f f(t) dt, 


T J-z/2 
7/2 
#22 f(t) cos ()a (re 1,2) 0) 
7 —1/2 vu 
7/2 
Boxe 1) sin (7 ar (n=1,2,...). 
T J-1/2 


(d) If desired, express the final Fourier series in a compact form with 
general formulas for its coefficients. 


To find the Fourier series for an odd periodic function f(t), do 
the following. 


(a) Identify f(t) as being odd, and find its fundamental period r. 
(b) Write down the Fourier series 
2° 
2nat 
F(t) =) Bnsin | —}. 
(t) > nsin ( = ) 


(c) Find the coefficients by evaluating the definite integrals 


4 7? _ (2nzt 
-if f(t) sin ( - at 


(d) If desired, express the final Fourier series in a compact form with 
general formulas for its coefficients. 
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To find the Fourier series for an even periodic function f(t), do 
the following. 


(a) Identify f(t) as being even, and find its fundamental period r. 


(b) Write down the Fourier series 
= 2nat 
F(t) = Ao +A cos (=) ‘ 
n= 


(c) Find the coefficients by evaluating the definite integrals 
2 7/2 
Ap == f(t) dt, 
T Jo 


7/ 
A, = <[ * Ft) eos (=) dt. 
T Jo 


T 


(d) If desired, express the final Fourier series in a compact form with 
general formulas for its coefficients. 


If f(t) is a periodic function, then any interval of length equal to its 
fundamental period 7 can be used to evaluate the Fourier coefficients. 


The following identities are often useful when evaluating Fourier 
series: if n is an integer, then 


cos(nm) = (—1)", 


sin(nz) = 0, 

(Mm) _ f (-1)"/? for n even, 
oon (F) = {6 for n odd, 
sin (=) _ fo for n even, 

2/7 | (-1)@*3)/2 for n odd. 


The following integrals may be useful for evaluating Fourier 
coefficients: 


[ tsin(at dt = + {sin(at) — at cos(at)] , 
J tcos(at) at = 4 {cos(at) + at sin(at)] . 


The even extension feyen(t) of a function f(t) defined on an interval 
[0,7] is the even periodic function with period 2T and fundamental 
interval [—T, T] defined by 


f(t) for0<t<T, 
Feven(8) = ine for -T <t <0, 
feven(t + 21) = feven(t)- 


14. 


16. 


Unit 12 summary 


The odd extension f,q(t) of a function f(t) defined on an interval 
[0, T] is the odd periodic function with period 2T and fundamental 
interval [—T,T] defined by 


= f(t) for OST, 
sane ee for -T <t <0, 


foaa(t + 2T) = foaa(t)- 


. If f(t) is a periodic function with Fourier series F(t), then F 


converges to f everywhere if f is continuous. If f(t) has a 
discontinuity at ¢ = to, then 


F(to) = 5lf (to) + f(t I, 
where f et ) is the limit of f(t) as t approaches tp from above, and 
f (to ) is the limit of f(t) as t approaches to from below. 


If a continuous periodic function f(t) with fundamental period 7 has 
Fourier series F(t), then its derivative f’(t) has the same fundamental 
period 7, and its Fourier series is given by F’(t). 


Unit 12 Partial differential equations 


Ls 


A partial differential equation is an equation relating a dependent 
variable and two or more independent variables through the partial 
derivatives of the dependent variable. The order of a partial 
differential equation is the order of the highest derivative that occurs 
in it. A partial differential equation is homogeneous if there are no 
terms that are solely functions of the independent variables. It is 
linear if the terms that contain the dependent variable are 
proportional to the dependent variable or to one of its partial 
derivatives. 


Examples of partial differential equations are the wave equation and 
the diffusion equation. 


The one-dimensional form of the wave equation is 
1 Pu Pu 


COP Ax’ 

where c is a positive constant called the wave speed. 

The one-dimensional diffusion equation is 

Ou 7 Ou 

Ot Ox? 

where D is a positive constant called the diffusion coefficient. 


The diffusion equation is also known as the heat equation. 
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Homogeneous linear partial differential equations satisfy the principle 
of superposition: if u(x,t) and u2(x,t) both satisfy a homogeneous 
linear partial differential equation, then any linear combination 


u(x,t) = a; uy (x, t) + a2 u2(x, t) 


(where a; and az are constants) also satisfies the partial differential 
equation. 


In order to find a unique solution of a partial differential equation for 
u(x,t), we must specify boundary conditions and initial conditions. 


A boundary condition is an equation relating to a specific value of 
a that holds for all time. For example, if an end of a taut string is 
held at a fixed position or a rod is held at a fixed temperature, then 
the boundary condition is u(L,t) = 0 for all t. If the end of a rod is 
thermally isolated, then the boundary condition is 


au 6,2) =0 for allt. 

Ox 
An initial condition is an equation that relates to a specific value 
of ¢ (usually ¢ = 0) and specifies u(x,t) at this time. For example, if 
the initial shape of a taut string (or the initial temperature profile of a 
rod) is given by f(x) for 0 < a < L, then the initial condition is 
u(x,0) = f(x) for 0 <a < L. If the initial velocity profile of a taut 
string is g(a) for 0 < x < L, then the initial condition is 

ee, 0) =g(r) fo0<2<L. 
To use the method of separation of variables to solve a linear 
homogeneous partial differential equation with dependent variable u 
and independent variables x and t, subject to given boundary and 
initial conditions, do the following. 


(a) Write the unknown function u(x,t) as a product of functions of 
one variable: 


u(x,t) = X(x) T(t). 


Find the required partial derivatives of u in terms of the ordinary 
derivatives of the functions X and T. For example, 
0u/da? = X"(x) T(t). 

(b) Substitute the partial derivatives found in step (a) into the partial 
differential equation. Rearrange the equation so that each side 
consists of a function of a single independent variable. Equate 
each side of the rearranged equation to the same separation 
constant j, and hence obtain ordinary differential equations for 
X and T. 


For example, the wave equation separates to give the equations 


X" (x) 1 T(t) 
xe)” era * 


(c) 


(d) 


(f) 


(g) 


Unit 12 summary 


Use the given boundary conditions for u to find boundary 
conditions for X. 


For example, the boundary condition u(0,t) = 0 for all ¢ will give 
the boundary condition X(0) = 0, and the boundary condition 
(du/Ox)(L,t) = 0 will give the boundary condition X’(L) = 0. 
Solve the differential equation for X, and apply the boundary 
conditions. Consider different choices for the separation 

constant yt. (Typically, the solutions X (x) take a different form 
depending on whether the separation constant is positive, 
negative or zero.) The boundary conditions generally produce a 
discrete set of solutions X,,(a), which are called eigenfunctions, 
and a corresponding discrete set of values j1,, for the separation 
constant, which are called eigenvalues. 


For each allowed j1,,, determine the corresponding solution T,,(#) 
of the differential equation for T(t). 


Combine X,,(x) and T;,(¢) to obtain a family of product solutions 
Un(x,t) = Xp(x) T(t), n=1,2,3,.... 


Express the general solution as an infinite linear combination of 
these product solutions: 


co 
u(x,t) = > an un(2,t). 
n=1 
Use the initial conditions and results about Fourier series to 


determine (when possible) the coefficients an. 


For example, the general solution for a taut string released from 
rest is 


+A ain (=) we nict 

ia L Li}? 
and if the initial profile of the string is given by the initial 
condition u(x,0) = f(x) for 0 < x < L, then the coefficients A, 


are given by the coefficients of the odd periodic extension of f(x), 
namely 


2). nar 
An == f sin (“*) f(v) ae, 10 12273 oc 


The general solution of the equation X”(x) = 4: X (zx) splits into three 
cases depending on whether ju is positive, negative or zero (all cases 
have two constants A and B). 


If 1 is positive, then js = c? for some c > 0, and the general 
solution is X(x) = Aexp(cr) + Bexp(—cr). 


If is negative, then yp = —k? for some k > 0, and the general 
solution is X(a) = Asin(kx) + Bcos(kx). 
If yz is zero, then the general solution is X(«) = Ar + B. 
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Unit 13. Non-linear differential equations 


1. A non-linear differential equation, or system of differential 
equations, is one that is not linear in at least one of the dependent 
variables. An autonomous differential equation, or system of 
differential equations, is one where the independent variable does not 
appear explicitly. 


This unit is mainly concerned with systems of two coupled 
autonomous differential equations, of the form 

dx dy 

—=u(x;y), = =v(z,y)- 

FH uey), B=) 

The system is non-linear if either u or v is a non-linear function. 


2. For systems of two coupled differential equations, the solution at 
time ¢t can be represented by a point x in a two-dimensional space 
called phase space, with coordinates «, y. 


The phase point x has equation of motion 
x=u(z,y), 
where u(x,y) is a vector field 


u=[u(z,y) v(x,y)]7. 
A particular solution [(t) y(t)]” of the equation of motion 
determines a path in the ry-plane, parametrised by t, whose tangent 
vector at any point (x,y) on it is the vector u(x,y). Such a solution 
curve is called a phase path (or phase trajectory). The xy-plane 
containing the solution curves is called the phase plane, and a 
diagram showing the phase paths is called a phase diagram. 

3. An equilibrium point (or fixed point) of a system of differential 
equations [¢ gy]? = u(z,y) is a point (ae, ye) such that x(t) = re, 
y(t) = ye is a constant solution of the system of differential equations; 
that is, ¢(t) = 0 and g(t) = 0 at the point (re, ye). 


To find the equilibrium points of x = u(x, y) for a given vector field u, 
solve the equation 

u(z,y) =0 
for x and y. 


4. An equilibrium point of a system of differential equations is said to be 
stable when all points in the neighbourhood of the equilibrium point 
remain in the neighbourhood of the equilibrium point as time 
increases, and unstable otherwise. 


5. The matrix 


Ou Ou 
an! y) ay y) 
J(z,y) = Ov Ov 


is the Jacobian matrix of the vector field u= [uv]. 
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Suppose that the system of differential equations x = u(x) has an 
equilibrium point at « = we, y = Ye. To linearise the system near this 


equilibrium point, do the following. 


(a) Find the Jacobian matrix J(a, y). 


(b) In the neighbourhood of the equilibrium point, the differential 


equations can be approximated by the linearised form 
Ou Ou 

| Fg (ter Ye) By Pere) Fy 
Ov Ov : 
Beltewe) Fp (ede) 


Unit 13 summary 


where x(t) = ze + p(t) and y(t) = ye + q(t). The vector [p q]7 is 


the perturbation from the equilibrium point. 


Consider a linear system of differential equations p= Ap. This 


system has a single equilibrium point, at the origin. The nature of the 
equilibrium point is determined by the eigenvalues and eigenvectors of 


the matrix A, as indicated by the following table. 


There are special cases not 
considered in this module (such 
as when an eigenvalue is zero). 


Real eigenvalues, 
distinct 


Both positive: source 


Both negative: sink 


Differing signs: saddle 


q 


Real eigenvalues, 


q 


equal 2 P 
Positive: star source Negative: star sink 
q 
q q 
Complex 
eigenvalues P P 


Positive real part: 


spiral source spiral sink 


Negative real part: 


Pure imaginary: 
centre 
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10. 


11. 


The equilibrium point is stable if it is a sink, a star sink, a spiral sink 
or a centre, and unstable if it is a source, a star source, a spiral 
source or a saddle point. 


To classify the equilibrium points of the non-linear system of 

differential equations x = u(x, y), do the following. 

(a) Find the equilibrium points using the method described in item 3. 

(b) Use the procedure described in item 6 to find the linear system 
that approximates the non-linear system in the neighbourhood of 
each equilibrium point. 

(c) Use the table in item 7 to classify the linear system, with the 
matrix A in item 7 equal to the linear approximation J(e, ye) 
found in item 6. 


The behaviour of a system of non-linear differential equations near an 
equilibrium point is the same as the behaviour of the linear 
approximation in the neighbourhood, except when the linear system 
has a centre. If the linear system has a centre, then the equilibrium 
point of the original non-linear system may be a stable centre, a stable 
spiral sink or an unstable spiral source. 
Some equations of motion have a constant of motion K (x,y). This 
function has the property that K (x(t), y(t)) remains constant along 
any given phase path (z(t), y(¢)). This implies that K satisfies a 
differential equation 
OK dx | OK dy _ 
Ox dt + Oy dt 
When there is a constant of motion, the phase paths must be contours 
of K(x,y). If the initial condition lies on a closed contour of K(x, y), 
then the motion will be periodic in t. 


0. 


An important application is population dynamics, where x(t) and 
y(t) are populations of two interacting species (usually predator and 
prey). The Lotka—Volterra equations are an important model: 


te ate(-2), Bo-m(1-3), 


where k, h, X and Y are positive constants. 


The equilibrium points of the Lotka—Volterra equations are an 
unstable saddle at (0,0) and a stable centre at (X,Y). 


The Lotka—Volterra equations have a constant of motion 


K(a,y) =hm«+kiny— te fy, 


The contours are all closed curves (except for the z- and y-axes). 


. A second-order differential equation # = f(x,<) can be converted into 


a pair of simultaneous first-order equations by setting y = #; the 
equivalent pair is 


t=y, y=f(x,y). 


13. 


14. 


Unit 13 summary 


The undamped pendulum satisfies the differential equations 


2 


f=y, y=-w'sing (-7<a2<7). 


This system has two equilibrium points: a stable centre at (0,0) and 
an unstable saddle at (7,0). 


The damped pendulum satisfies the differential equations 
é=y, y= -w'sing—ey (—1<2<7n), 


where ¢ is a positive constant. This system has two equilibrium points, 
at (0,0) and (7,0). The equilibrium point (0,0) is a stable spiral sink 
for 0 < ¢ < 2w, a stable improper sink for ¢ = 2w, and a stable sink for 
€ > 2w. The equilibrium point (7,0) is an unstable saddle. 
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Index 


acceleration 19 
acute angle 13 
addition 
of matrices 36 
of vectors 33 
amplitude of a sinusoidal function 15, 30 
analytic methods of solution 28 
angle 
between vectors 34 
subtended by an arc of acircle 18 
angle sum 
of a quadrilateral 17 
of atriangle 17 
angular frequency of a sinusoidal function 15, 30, 64 
antiderivative 24 
antiparallel vectors 33, 35 
arbitrary constant 
of a differential equation 27 
of integration 24 
arc length 18 
arc of acircle 18 
arccos 14 
arccosec 14 
arccosh 11 


arccot 14 
aresec 14 
aresin 14 
aresinh = 11 
arctan 14 
arctanh 11 
area 


ofacircle 18 
of a parallelogram 17, 35 
of asector 18 
ofatriangle 17, 35 
under a graph 24 
area element 
in orthogonal coordinates 56 
in polar coordinates 54 
on a curved surface 57 
area integral 
and calculation of area 54 
in polar coordinates 54 
of a product function with constant limits 52 
over arectangle 52 
over an arbitrary region in Cartesian coordinates 
53 
Argand diagram 8 
argument of a complex number 8 
associated homogeneous equation 31 
associativity of matrix addition 36 
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asymptote 23 
augmented matrix 41 


autonomous differential equation 70 


auxiliary equation 30 
axially-symmetric region 


back substitution 42 


55 


base of an exponential function 10 


basis 43 


boundary condition 32, 


boundary value 32 
boundary-value problem 
bounded above 22 
bounded below 22 


68 


32 


Cartesian components of a vector 


Cartesian coordinates 


33. 


for three-dimensional space 33 


for two-dimensional 


space 16 


Cartesian form of a complex number 


Cartesian unit vector 3. 
centre of a circle 18 


3, 40 


chain rule of differentiation 20 
chain rule of partial differentiation 
for a change of variables 49 


for differentiation with respect to a parameter 


for small changes 49 
characteristic equation of a matrix 44 


chord 
ofacircle 18 
ofacurve 18 
circle 18 
circulation 63 
circumference of a circle 


18 


classification of a stationary point 


using eigenvalues of the Hessian matrix 51 


using the determinant test 51 
classification of an equilibrium point 


closed interval 7 
closed surface 62 

unit normals of 62 
coefficient matrix 41 


coincident roots of a polynomial 9 


column vector 36 


commutativity of matrix addition 36 


complementary function 


31 


7 


71, 72 


of a system of differential equations 47 


complex conjugate 7 
complex number 7 
complex-valued function 


19 


component form of a vector 33-35 


component of a vector 


33, 34 


49 


composite function 9 
composite rule for differentiation 20 
composition of functions 9 
conservative field 61 
constant 8 
of integration 24 
constant function 9 
constant multiple rule 
for differentiation 20 
for integration 25 
constant of motion 72 
continuous function 9 
contour line 48 
contour map 48 
contour surface 48 
convergence 
of a sequence of numbers 6 
of Fourier series 67 
coordinate line 55 
corresponding homogeneous system 47 
cos 13 


cosec 13 
cosh 11 
cot 13 


cross product of vectors 35 

cubic function 9 

cuboid 52 

curl 
as circulation per unit area 63 
in Cartesian coordinates 59 
in cylindrical coordinates 59 
in polar coordinates 59 
in spherical coordinates 59 
interpretation of 59 

curl test 62 

curl theorem 63 

curve sketching 22 

curved surface 
area element 57 
area of 57 

cyclic identity 39 

cylindrical coordinates 54 


damped pendulum 73 
de Moivre’s theorem 8 
decimal places 6 
decreasing function 22 
definite integral 24 
degree of a polynomial 9 
degrees as a measure of angle 13 
dependent variable of a function 8 
derivative 19 
of a complex-valued function 19 
derived function 19 
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determinant 38 
ofa 2x2 matrix 38 
ofa3x3 matrix 38 
of ann x n matrix 38 
properties of 39 
rule for eigenvalues 45 
determinant test 51 
diagonal matrix 41 
diameter of a circle 18 
difference of two squares 10 
differential equation 27 
differentiation 19 
tules 20 
diffusion coefficient 67 
diffusion equation 67 
dilation matrix 40 
dimension 43 
direct integration 28 
direction of a vector 32 
discriminant of a quadratic equation 9 
displacement 32 
vector 32, 34 
distinct eigenvalues 43 
distinct eigenvectors 43 
distributivity of matrix addition and scaling 36 
divergence 
as flux per unit volume 63 
in Cartesian coordinates 59 
in cylindrical coordinates 59 
in polar coordinates 59 
in spherical coordinates 59 
interpretation of 59 
divergence theorem 63 
dot product of vectors 34 


eigenfunction 69 
eigenvalue 43, 69 
of a symmetric matrix 45 
of a triangular matrix 45 
of an invertible matrix 45 
eigenvector 43 
equations 44 
expansion 43 
element of a matrix 35 
ellipse 18 
equal roots of a polynomial 9 
equality 
of matrices 36 
of vectors 32 
equation of a circle 18 
equilateral triangle 13, 17 
equilibrium point of a system of differential equations 
70, 71 
Euler's formula 8 
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even extension of a function 66 

even function 64 

explicit solution of a differential equation 27 
exponent of an exponential function 10 
exponential behaviour 10 

exponential form of a complex number 8 
exponential function 10 


factorial 23 
factorisation of polynomials 10 
field 57 
first-order ordinary differential equation 27 
first-order partial derivative 49 
first-order system of differential equations 46 
first-order Taylor polynomial 23, 50 
fixed point 70 
flux of a vector field 

over a planar element 62 

over an extended surface 62 
formula method for roots of a quadratic equation 9 
Fourier coefficients 65 
Fourier series 65 

for a periodic function 65, 66 

for an even periodic function 66 

for an odd periodic function 65 
function 8 
function notation (for derivatives) 19 
function of a function rule for differentiation 20 
fundamental interval of a periodic function 64 
fundamental period of a function 64 


Gaussian elimination 42 
Gauss’s theorem 63 
general solution 
first-order, homogeneous 46 
first-order, inhomogeneous 47 
of a differential equation 27 
of a first-order system of differential equations 46 
of a homogeneous linear constant-coefficient 
second-order differential equation 30 
of a homogeneous second-order system of differential 
equations 48 
of an inhomogeneous linear differential equation 31 
global maximum of a function 22 
global minimum of a function 22 
gradient 58 
in Cartesian coordinates 58 
in cylindrical coordinates 58 
in polar coordinates 58 
in spherical coordinates 58 
properties of 58 
gradient (slope) 22 
gradient field 61 
gradient vector 50 
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and contour lines 50 
and slope 50 
gradient vector field 58 
graph of a function 9 
graphs of common functions 11-12 
Greek alphabet 5 


half-open interval 7 
heat equation 67 
Hessian matrix 51 
hexagon 17 
higher derivatives 19 
higher-order partial derivative 49 
homogeneous linear differential equation 
first-order 28 
second-order 29 
homogeneous partial differential equation 67 
homogeneous system 
of first-order differential equations 46 
of second-order differential equations 47 
horizontal axis of a Cartesian coordinate system 
hyperbola 18 
hyperbolic functions 11 
hypotenuse 14 


identity matrix 37 
imaginary part of a complex number 7 
implicit differentiation 20 
implicit solution of a differential equation 27 
inconsistent system 42 
increasing function 22 
indefinite integral 24 
independent variable of a function 8 
index of an exponential function 10 
inhomogeneous linear differential equation 
first-order 28 
second-order 29 
inhomogeneous system of first-order differential 
equations 46 
initial condition 68 
for a first-order differential equation 27 
for a second-order differential equation 32 
initial value 
for a first-order differential equation 27 
for a second-order differential equation 32 
initial-value problem 
for a first-order differential equation 27 
for a second-order differential equation 32 
integer 6 
integral 24 
integrand 25 
integrating factor 29 


method 29 
integration 25 
by parts 25 
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by substitution 25 

tules 25 
interval 7 
inverse function 10 
inverse hyperbolic functions 11 
inverse of a matrix 38 

2x2 formula 38 

nxn procedure 39 
inverse trigonometric functions 14 
invertible matrix 38 
irrational number 6 
isosceles triangle 17 


Jacobian matrix 70 
Jacobian vector 56 


Laplace’s rule 38 
leading diagonal 41 
Leibniz notation (for derivatives) 19 
length of a path 60 
limit of a sequence of numbers 6 
line integral 
of a gradient field 61 
of a scalar function 60 
of a vector field 60 
path-independence 61 
linear dependence 42 
linear differential equation 
constant-coefficient 29 
first-order 28 
second-order 29 
linear equation 41 
linear function 9 
linear independence 43 
linear partial differential equation 67 
linear transformation 40 
linearised form of differential equations 
local maximum 22, 51 
local minimum 22, 50 
logarithm function 10 
Lotka-Volterra equations 72 
lower bound of a function 22 
lower limit of integration 24 
lower triangular matrix 41 


magnitude of a vector 32 
main diagonal 41 


matrix 35 
addition 36 
element 35 
inverse 38 


inversion 39 
multiplication 37 
of coefficients 41, 46 
order 35 
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maximum of a function 22 

method of undetermined coefficients 31 
minimum of a function 22 

mixed partial derivative theorem 49 
modulus of a complex number 7 


natural logarithm function 10 

Newtonian notation (for time derivatives) 19 
Newton's notation (for time derivatives) 19 
n-gon 17 

non-invertible matrix 38 

non-linear differential equation 70 
non-singular matrix 38 

normal mode of oscillation 48 

nth derivative 19 

nth-order polynomial 9 

nth-order Taylor approximation 23 
nth-order Taylor polynomial 23 


obtuse angle 13 
odd extension of a function 67 
odd function 25, 64 
open interval 7 
open surface 62 
order 
of a derivative 19 
of a differential equation 27 
of a matrix 35 
of a partial differential equation 67 
ofa polynomial 9 
of a Taylor approximation 23 
of a Taylor polynomial 23 
oriented area 62 
origin of a Cartesian coordinate system 16 
orthogonal coordinate system 55 
oscillatory function 16 


parabola 9, 18 
parallel vectors 35 
parallelepiped 35, 39 
parallelogram 17, 35 
parametric representation of a path 60 
partial derivative 48 
partial differential equation 67 
particular integral 
of a system of differential equations 47 
of an inhomogeneous differential equation 31 
particular solution 
of a differential equation 27 
of a first-order system of differential equations 46 
path 60 
parametric representation of 60 
path-independence of line integrals 61 
pendulum equation 73 
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pentagon 17 
perfect square 10 
period 
of a function 64 
of a sinusoidal function 16, 30 
periodic function 16, 64 
perpendicular vectors 34 
perturbation from an equilibrium point 71 
phase constant of a sinusoidal function 15, 30 
phase diagram 70 
phase path 70 
phase plane 70 
phase point 70 
phase space 70 
phase trajectory 70 
point of inflection 22 
polar coordinates 17, 54 
polygon 17 
polynomial 9 
function 9 
of degreen 9 
population dynamics 72 
position vector 34 
power 
of a matrix 37 
of an exponential function 10 
principal value 
of a polar coordinate angle 17 
of the argument of a complex number 8 
principle of superposition 30, 68 
product form 52 
product of matrices 37 
product rule for differentiation 20 
Pythagoras’s theorem 14 


quadrants of the plane 16 
quadratic approximation 23 
quadratic function 9 

quadrilateral 17 

quotient rule for differentiation 20 


radian 13 

radius of a circle 18 

rational number 6 

real matrix 45 

real number 6 

real part of a complex number 7 
rectangle 17 

regular polygon 17 

repeated eigenvalues 43 
repeated root of a polynomial 9 
resolving a vector 33 

resultant of vectors 33 

right angle 13 
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right-angled triangle 17 

right-hand grip rule 63 

right-hand rule for a coordinate system 33 
right-hand-side vector 41 

right-handed coordinate system 33 

root of a polynomial 9 

rotation matrix 40 

rounding 7 

row vector 36 


saddle point 51 
scalar 32 
scalar field 57 
scalar multiple 32 
of a matrix 36 
scalar multiplication 32 
scalar potential field 61 
scalar product of vectors 34 
scalar triple product of vectors 39 
scale factor 56 
scaling a vector 32 
scientific notation 7 
sec 13 
second derivative 19 
second-order ordinary differential equation 29 
second-order partial derivative 49 
second-order system of differential equations 47 
second-order Taylor polynomial 23, 50 
section function 48 
separation constant 68 
separation of variables 28, 68 
significant figures 6 
simultaneous linear equations 41 
sin 13 
singular equations 42 
singular matrix 38 
singular system 42 
sinh 11 
sinusoid 15 
sinusoidal function 15 
slope 
of a straight line 9 
of asurface 50 
small-angle approximations 24 
smooth function 22 
solution 
of a differential equation 27 
of a homogeneous differential equation 30 
of an inhomogeneous differential equation 31 
spherical coordinates 55 
square 17 
square matrix 36 
stable equilibrium point 70, 72 
standard derivatives 21 


standard integrals 25 
stationary point 50 
classification using the determinant test 51 
classification using the Hessian matrix 51 
local maximum 51 
local minimum 50 
local saddle point 51 
of a function of one variable 22 
Stokes’s theorem 63 
straight line, vector equation of 34 
subtraction of vectors 33 
sum 
of functions 8 
of matrices 36 
of vectors 33 
sum rule 
for differentiation 20 
for integration 25 
summation notation 6 
surface 48 
symmetric matrix 45 
system 
of differential equations 46 
of linear equations 41 
of linear first-order differential equations 46 
of linear second-order differential equations 47 
of non-linear differential equations 70, 72 


tan 13 
tangent approximation 23 
tangent plane 50 
tangent to a curve 18 
tanh 11 
Taylor approximation 23 
Taylor polynomial 
first-order 50 
of a function of one variable 23 
second-order 50 
Taylor series of a function of one variable 23 
third derivative 19 
trace of a matrix 44 
rule for eigenvalues 45 
transformation matrix 40 
transpose of a matrix 37 
trial solution 
for a linear differential equation 31 
for a system of differential equations 47 
triangle 17, 35 
triangle rule for vector addition 33 
triangular matrix 41 
trigonometric addition formulas 15 
trigonometric double-angle formulas 15 
trigonometric functions 13 
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trigonometric identities 14 


undamped pendulum 73 
uniqueness of solutions of a first-order differential 
equation 28 
unit vectors 33, 57 
in Cartesian coordinates 57 
in cylindrical coordinates 57 
in polar coordinates 57 
in spherical coordinates 57 
unstable equilibrium point 70 
upper bound of a function 22 
upper limit of integration 24 
upper triangular matrix 41 


variable 8 
vector 32 
vector addition 33 
vector equation of a straight line 34 
vector field 57 
conversions 58 
vector product of vectors 35, 39 
vector space 43 
vector subtraction 33 
velocity 19 
vertical axis of a Cartesian coordinate system 16 
volume 
of a parallelepiped 35, 39 
of an axially-symmetric region 55 
volume element 
in cylindrical coordinates 54 
in orthogonal coordinates 56 
in spherical coordinates 55 
volume integral 
and calculation of volume 54 
in cylindrical coordinates 55 
in spherical coordinates 55 
of a product function with constant limits 52 
over a cuboid 52 
over an arbitrary region in Cartesian coordinates 
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wave equation 67 
wave speed 67 


a-axis of a Cartesian coordinate system 16 
ay-plane 16 


y-axis of a Cartesian coordinate system 16 
y-intercept of a straight line 9 


zero function 9 
zero matrix 36 
zero vector 32 
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